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NOTE 


The  material  in  this  Report  is  intended  to  supplement  and  provide 
a more  detailed  fluid  dynamic  input  to  AGARD  Advisory  Report 
No.  47,  The  Physical  Vulnerahilitv  of  Aircraft,  (3  volumes,  Sept- 
ember 1972,  May  1973  and  May  1973,  NATO  Classified). 
Readers  with  a general  interest  in  the  subject  of  aircraft  physical 
vulnerability  are  referred  to  AR-47. 
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This  Advisory  Report  has  been  prepared  at  the  request  of  the  Fluid  Dynamics  Panel  of  AG.ARD. 


THE  MISSION  OF  AGARD 


The  mission  of  AGARD  is  to  bring  together  the  leading  personalities  of  the  NATO  nations  in  the  fields  of 
science  and  technology  relating  to  aerospace  for  the  following  purposes: 

- Exchanging  of  scientific  and  technical  information; 

- Continuously  stimulating  advances  in  the  aerospace  sciences  relevant  to  strengthening  the  common  defence 
posture; 

- Improving  the  co-operation  among  member  nations  in  aerospace  research  and  development; 


- Providing  scientific  and  technical  advice  and  assistance  to  the  North  Atlantic  Military  Committee  in  the 
field  of  aerospace  research  and  development; 

- Rendering  scientific  and  technical  assistance,  as  requested,  to  other  NATO  bodies  and  to  member  nations 
in  connection  with  research  and  development  problems  in  the  aerospace  field; 

- Providing  assistance  to  member  nations  for  the  purpose  of  increasing  their  scientific  and  technical  potential; 

- Recommending  effective  ways  for  the  member  nations  to  use  their  research  and  development  capabilities 
for  the  common  benefit  of  the  NATO  community. 


The  highest  authority  within  AGARD  is  the  National  Delegates  Board  consisting  of  officially  appointed  senior 
representatives  from  each  member  nation.  The  mission  of  AGARD  is  carried  out  through  the  Panels  which  are 
composed  of  experts  appointed  by  the  National  Delegates,  the  Consultant  and  Exchange  Program  and  the  ,\erospace 
Applications  Studies  Program.  The  results  of  AGARD  work  are  reported  to  the  member  nations  and  the  NATO 
Authorities  through  the  AGARD  series  of  publications  of  which  this  is  one. 
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TtCHNIQUtS  FOR  PRtOICTING  STRUCTURAL  RI  SPONSt 
DUt  TO  tXPLOSIVE  AIR  BLAST 


INTRODUCTION 


The  problem  of  predicting  the  response  of  an  aircraft  structure  to  external  blast  Is  difficult  for 
the  simple  reason  that  the  aircraft  designer  has  never  been  required  to  design  aircraft  that  would 
survive  the  blast  wave  of  a warhead.  Since  the  aircraft  designer  has  never  been  requested  to  design 
aircraft  that  could  survive  a blast  wave,  the  aircraft  stress  analyst  has  never  had  to  develop  the 
tools  for  predicting  the  response  of  aircraft  skin  and  structure  to  a free  air  blast. 

On  the  other  hand,  the  vulnerability  analyst  and  the  warhead  designer  are  continually  consulted 
about  the  size  of  warhead  needed  to  defeat  enemy  aircraft.  They  are  prevented  from  taking  a ratlr.ial 
approach  and  really  settling  the  problem  by  at  least  three  factors: 

1.  Construction  details  of  new  enemy  aircraft  are  obviously  unavailable  to  them.  Usuall),  only 
some  photographs  of  the  external  configuration  are  available.  Hence,  a rational  structural  response 
model  from  first  principles  (i.e.,  Newton's  laws,  etc.)  cannot  be  developed  since  they  do  not  have 
access  to  the  necessary  details  to  complete  a rational  model. 

2.  The  main  kill  mechanism  of  warheads  is  usually  fragmentation,  with  air  blast  being  a secondary 
kill  mechanism.  Hence,  most  work  has  been  done  on  fragmentation  effects  modeling. 


3.  The  vulnerability  analyst  and  the  warhead  designer  are  usually  quite 
"nuts  and  bolts"  approach  of  the  aircraft  stress  analyst  and,  therefore,  they 
engineering  language  that  the  aircraft  designer  and  stress  analyst  use.  This 
fact  that  the  aircraft  designer  has  never  had  to  consider  warhead  effects  and 
warheads  and  structural  response  to  blast  loading. 
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APPROACH 

The  approach  that  will  be  taken  to  the  problem  of  predicting  the  structural  response  to  external 
blast  is  that  of  adapting  and  summarizing  the  most  important  details  of  a spherical  charge  in  free  air 
as  it  Interacts  with  rigid  surfaces. 

The  use  of  impulse  momentum  foi  calculating  structural  response  and  the  approx imat ion  of  a uni- 
formly distributed  impulse  on  tlie  surface  in  question  are  approximations  that  appear  to  be  experimen- 
tally justified. 

Also  listed  is  a technique  by  Sewell  and  Kinney  (Crjtjcai  J/npuise  K’lthi/j  a critical  Time),  which 
has  been  used  by  designers. 


SAFETY  FACTOR  AND  ACCURACY 

Certain  simplifications  have  been  made  in  the  development  of  the  design  procedures  presented  in 
this  manual.  As  a result,  an  analysis  of  a structure  using  these  procedures  will  generally  result  in 
a conservative  estimate  of  the  structure's  capacity  and,  consequently,  structures  designed  using  these 
procedures  will  generally  be  adequate  for  the  blast  load  exceeding  the  assumed  loading  conditions. 

Certain  unknown  factors  can  result  in  an  over-estimate  of  the  protective  structure's  capability 
to  resist  the  effects  of  an  explosion.  These  factors,  reflections  of  the  shock  waves,  structural 
response,  construction  methods,  quality  of  construction  materials,  workmanship,  etc.,  vary  for  each 
design.  To  compensate  for  weaknesses  resulting  from  these  factors,  it  is  recommended  that  the  "effec- 
tive charge  weight"  or  the  actual  charge  weight,  depending  upon  the  method  used  to  determine  the  TNT 
equivalent,  be  increased  by  20Z  for  design  purposes. 


The  effects  of  the  warhead  shape  (which  can  increase  the  Impulse  and  peak  pressure)  are  neglected. 
The  effects  created  when  the  warhead  is  in  direct  contact  with  the  structure  are  also  neglected 
(neglected  here  because  at  the  present  time,  they  are  Impossible  to  calculate). 

Another  directional  effect  is  sometimes  found  in  high-velocity  munitions  (such  as  23-mm  high- 
explosive  incendiary  (UEl)  , etc.),  particularly  in  contact  or  close-in  proximity  fuzing.  Here,  the 
forward  velocity  of  the  munitions  will  be  felt  in  an  increased  damage  capability  in  the  forward  direc- 
tion. Use  of  impulse  momentum  structural  response  technique  infers  that  the  momentum  of  fragment  and 
gas  cloud  be  added  to  the  blast  wave  Impulse.  This  effect,  to  the  writer's  knowledge,  has  only  been 
observed  qualitati''ely , and  has  received  little  quantitative  verification. 

Another  effect  which  has  received  little  attention  is  that  of  warhead-aircraft  closing  velocity. 
This  has  been  discussed  at  length,  but  no  hard  and  fast  rules  as  to  how  to  include  it  in  warhead  air 
blast  effects  have  been  formulated. 


COMMENTS  ON  DATA  SOURCE 

The  data  source  listed  in  Reference  1 was  selected  because  the  approach  taken  was  a rational  one 
and  the  theoretical  predictions  were  experimentally  verified.  Although  the  results  quoted  in  the 
reference  were  obtained  with  concrete  walls,  the  writer  used  the  same  technique  on  walls  built  of  steel 
and  obtained  good  agreement  between  theory  and  experiment.  In  terms  of  pure  blast  loading  of  struc- 
tures, the  civil  engineering  profession  has  had  the  most  demands  for  rational  and  useful  design  proce- 
dures for  calculating  the  effects  of  air  blast  loading  on  structures.  The  advent  of  nuclear  weapons 
has  forced  the  civil  engineer  to  deal  with  the  design  of  a building  to  resist  air  blast.  The  problems 
of  designing  explosive  storage  and  loading  facilities  so  that  catastrophic  incidents  would  not  occur 
has  also  forced  the  civil  engineer  to  be  able  to  rationally  design  structural  walls  that  can  withstand 
blast  loads.  The  closest  approach  in  tlie  field  of  damage  to  aircraft  structures  has  been  tiie  "critical 
Impulse  in  a critical  time"  technique.  Other  characterization  such  as  the  generic  kill  damage  contours 
is  of  little  value  to  the  aircraft  designer  since  it  does  not  relate  improvements  in  design  in  their 
ability  to  resist  attack. 


OUTLINE  OF  FUNDAMENTALS 

1.  Spherical  charges  in  air,  Hopklnson  and  Sachs'  scaling  relations 

2.  Reflected  pressure  coefficient  versus  angle  of  incidence 

3.  External  blast  loads  on  structures 

a.  Reflected  pressure 

b.  Clearing  time 

c.  Drag  phase 

d.  Geometric  effects 

4.  Impulse  momentum  approach 

5.  Critical  Impulse  in  a critical  time 

6.  Closed  compartment  estimates 

7.  Summary  and  limitations 


SPHERICAL  CHARGES  IN  AIR 


Blast  Wave  Phenomena 

The  violent  release  of  energy  from  a detonation  in  a gaseous  medium  gives  rise  to  a sudden  pressure 
Increase  In  that  medium.  The  pressure  disturbance,  termed  the  blast  wave,  Is  characterized  by  an  almost 
instantaneous  rise  from  the  ambient  pressure  to  a peak  incident  pressure  ^so‘  This  pressure  increase 
or  shock  front  travels  radially  from  the  burst  point  with  a diminishing  velocity  U,  which  is  always  In 
excess  of  the  sonic  velocity  of  the  medium.  Gas  molecules,  nwiking  up  the  front,  move  at  lower  veloc- 
ities u.  This  latter  particle  velocity  is  associated  with  a dynamic  pressure  or  the  pressure  formed 
by  the  winds  produced  by  the  shock  fronts.  As  the  shock  front  expands  into  increasingly  larger  volumes 
of  the  medium,  the  peak  incident  pressure  at  the  front  decreases  and  the  duration  of  the  pressures 
Increases . 


3 


At  any  point  away  from  the  burst,  the  pressure  disturbance  has  the  shape  shown  in  Figure  1.  The 
shock  front  arrives  at  time  and,  after  the  rise  to  the  peak  value,  the  incident  pressure  decays  to 
the  ambient  value  in  the  time  which  is  the  positive  phase  duration.  This  is  followed  by  a negative 
phase  with  a duration  t©”  longer  than  the  positive  phase  and  characterized  by  a pressure  below  the 
pre-shot  ambient  pressure  (maximum  value  of  Pso~^  ^ reversal  of  the  particle  flow.  The  negative 
phase  is  usually  less  important  in  a design  than  is  the  positive  phase.  The  Incident  impulse  density 
(hereafter  referred  to  as  unit  Incident  Impulse)  associated  with  the  blast  wave  is  the  integrated 
area  under  the  pressure-time  curve  and  is  denoted  as  for  the  positive  phase  and  ig”  for  tl>e  negative 
phase . 


Fig.)  Free  field  pressure  tune  variation 


If  the  shock  wave  impinges  on  a rigid  surface  oriented  at  an  angle  to  the  direction  of  propagation 
of  the  wave,  a reflected  pressure  is  instantly  developed  on  the  surface,  and  the  pressure  Is  raised 
to  a value  In  excess  of  the  incident  pressure.  The  reflected  pressure  is  a function  of  the  pressure 
in  the  incident  wave  and  the  angle  formed  between  the  rigid  surface  and  the  plane  of  the  shock  front. 

For  a reflector,  where  flow  around  an  edge  or  edges  occurs,  the  duration  of  the  reflected  pressures 
is  controlled  by  the  size  of  the  reflecting  surface.  The  high  reflected  pressure  seeks  relief  toward 
the  lower  pressure  regions,  and  this  tendency  is  satisfied  by  the  propagation  of  rarefaction  waves 
from  the  low-  to  the  high-pressure  region.  These  waves,  traveling  at  the  velocity  of  sound  in  the 
reflected  pressure  region,  reduce  the  reflected  pressures  to  the  stagnation  pressure  which  is  the  value 
that  is  in  equilibrium  with  the  high-velocity  air  stream  associated  with  the  incident  pressure  wave. 

When  such  a relief  is  not  possible  (for  example,  when  an  incident  wave  strikes  an  infinite  surface), 
the  Incident  pressure  at  every  point  in  the  wave  will  be  reflected,  and  these  reflected  pressures  will 
last  for  the  duration  of  the  wave. 

The  peak  positive  reflected  pressure  is  designated  P^-,  the  peak  negative  reflected  pressure  Pj*“, 
and  the  unit  impulses  associated  with  a completely  reflected  incident  wave  are  for  the  positive 
phase  and  i for  the  negative  phase  (Figure  2 lists  the  behavior  of  a spherical  charge  at  sea  level.) 

Departures  From  Spherical  Charge  Behavior 

For  warheads  other  than  spherical  shapes,  the  blast  wave  produced  may  not  be  spherically  synanetri- 
cal.  Figure  3 (from  Reference  2)  illustrates  the  variation  of  the  peak  pressure  (and  impulse)  of  a 
point- initiated  cylindrical  charge.  The  pressure  and  impulse  may  be  influenced  to  a much  greater  extent 
by  the  use  of  multiple  detonators,  focusing  Jacket,  etc.  It  is  recommended  that  the  designer  very 
carefully  inquire  as  to  the  extent  of  the  nonspherical  behavior  of  the  threat  warhead  before  he  attempts 
a blast-induced  structural  response  calculation. 
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’ ^AK  fOSiriVE  IMCIOCMT  MESSuAf  m 
• PIAK  NEGATIVE  INCIDENT  f«CSSURl  m 
f,  • AEAK  AOSiriVE  NONMAl  NEFLECTEO  FNESSUKE 
f;  • FEAK  NEGATIVE  NOAMAl  NEFlECTEO  MESSUftf  ^ 

• SCACEO  UNIT  fDtiTiVE  INCIDENT  iMtuiSE  P%\  W IS'  ^ 

• ICAIEO  UNIT  NEGATIVE  INCIDENT  impulse  PSi  MS  IB' 3 

• SCALED  UNIT  POSITIVE  NORMAL  REFLECTED  IMPULSE  PSi  MS  LB*  3 

• SCALED  UNIT  NEGATIVE  NONMAL  REFLECTED  IMPULSE  PS*  MS  LB'  3 
ia/Mi<'3  - SCALED  TIME  OF  ARRIVAL  OF  BLAST  MAVE  Ml.ll*^ 


ioW<  3 • scaled  positive  duration  OF  POSiTivE  PHASE  MS  IB*'3 
■ oTR'  3 ■ SCALED  NEGATIVE  DURATION  OF  POSITIVE  PHASE  MS  LB’  3 
L«  W’  3 • SCALED  AAVE  LENGTH  OF  POSiTivE  PHASE  FT  iB>  3 
l;;W'  3 < SCALED  WAVE  LFNGTH  OF  NEGATIVE  PHASE  FT  L#'  3 
U • SHOCK  front  velocity  FT  ms 
u - PARTICLE  VELOCITY  FT  MS 
W charge  WEIGHT  LB 

R RACIAL  Distance  FROM  LHAROE  ft 
/ • scaled  Distance  ft  ib<  3 
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Fjg.3  Ratio  of  free  air  peak  overpressure  (P-cylindcr/P-sphcre)  versus  distance  for  cylinders  with  differing  aspect  ratios  (i'/d) 


TNT  Kqui valent.  The  major  quantity  of  blast  effects  data  presented  in  this  manual  pertains  to  the 
blast  pressure  output  of  TNT  explosions.  These  data  can  be  extended  to  Include  other  potentially  mass- 
detonating  materials  and  explosives  whose  shapes  differ  from  those  considered  by  relating  the  explosive 
energy  of  the  effective  charge  weight  of  these  materials  to  that  of  an  equivalent  weight  of  TNT.  The 
explosive  energy  of  the  effective  charge  weight  is  also  equal  to  that  of  the  total  mass  of  the 
material.  No  clear-cut,  high-explosive  equivalent  may  be  utilized  to  define  the  blast  effects  of 
explosive  materials  other  than  high  explosives  themselves.  The  blast  effects  of  other  materials  in 
the  anticipated  environments  must  be  analyzed  and  then  related  to  the  blast  effects  produced  by  the 
TNT  explosion  at  the  range  of  interest  to  obtain  their  equivalencies.  To  illustrate  a typical  analysis, 
explosive  tests  of  certain  propellant  liquids  and  hydrocarbon  mixtures  have  indicated  that  the  explo- 
sive equivalent  of  these  materials,  which  relate  both  the  peak  blast  pressures  and  impulse,  is  constant 
over  the  entire  intermediate-  and  low-pressure  ranges.  At  higher  pressures,  the  TNT  equivalent  will 
vary  for  each  pressure  level  and  will  be  different  than  the  TNT  equivalent  which  relates  the  impulse. 


For  blast-resistant  design  in  general,  the  TNT  equivalent  should  be  based  upon  a pressure  and/or 
impulse  relationship  depending  upon  the  anticipated  pressure-design  range.  In  addition  to  being  a 
function  of  the  comparable  blast  effects  themselves,  the  equivalent  will  also  be  dependent  upon  the 
anticipated  structural  response.  For  the  lower  pressure  ranges  where  structures  respond  to  either  the 
peak  pressure  or  pressure  pulse  of  the  shock  wave,  the  ratio  of  the  weights  of  given  explosives  to 
that  of  TNT  is  given  in  Table  1 for  both  peak  pressure  and  impulse  in  pressure  levels  varying  between 
2 and  50  psi. 


LONGITUDINAL  AXIS 
CYLINDER 
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Comp  A- 3 
Comp  B 

Comp  B/riWT  (70/30) 

Cyclotol  (70/30) 

Explosive  D 

HBX-1 

HBX-3 

H-6 

Minol  II 


tabu;  1 

tquivalem  Weight  Ratios  for  Free-Air  Effects 


Pentolite 
Picratol 
RDX/Wax  (98/2) 
RDX/Wax  (95/5) 
TNT 
TNETB 
Torpex  II 
Tritonal  (80/20) 


Peak  pressure 
Per, 


NOTE:  Data  are  applicable  for  pressure  levels  between  2-  and  50*-psi  range  for  shock  overpressure 

HOPKINSON  AND  SACHS’  SCALING  LAWS 

The  data  in  Figure  2 are  scaled  for  sea  level  conditions  according  to  Hopkinson's  law.  Another 
scaling  law  is  Sachs'  law.  Figures  4,  5,  and  6 (from  Reference  3)  show  peak  pressure  and  scaled 
impulse  for  a range  of  conditions.  Of  interest  is  the  size  of  experimental  error. 
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hg.4  Peak  overpressure  versus  Sachs-scaled  distance  (Ref  4) 
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Fig. 5 Hopkinson-scaleU  impulse  versus  distance  (Ref.  4) 
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standard  deviation  of  aeasure- 
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( ig.(>  Sachs  scaled  impulse  versus  distance  (Ret  41 


Sachs'  1 .iw  lias  boon  conflrmod  by  oxporimcnis  «*  Uowt-y  and  Sporra;tza,‘  Iricsson  and  Kdin,  and 
of  Olson,  ot  al.,*  with  by  tar  tho  most  oxtonsivo  scrlos  ot  iiK’dt'l-prolotypc  ct'inparisons  being  Lho 
firsL.  Dewoy  .ind  Sporrazza**  conJin  tod  thoir  tests  with  several  sizes  of  bare  Pentolile  splieros  in 
an  a 1 1 1 1 ude-s i mu  1 .U  i ng  chamber  in  which  both  the  .unbient  pressure  and  temperature  could  be  varied. 
Arrays  ot  side-on  blast  pressure  transducers  were  mounted  at  various  distances  from  the  expli'sive 
spheres,  and  time  histi’ries  ot  pr»-ssurf  were  recorded.  The  two  priniary  blast  pariuneters  re/orted 
were  peak  overpressure  f .ind  positive  impulse  I.  A number  ot  repeal  tests  were  ci'iidui  ted  lor  any 
given  condition  and  distance.  Figure  U shows  their  dal.i  tor  peak  overpn-ssure , scaled  according  t 
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Sachs'  law.  Their  data  for  impulse  with  Hopkinson  scaling  applied  are  shown  in  Figure  5,  and  with 
Sachs'  scaling,  in  Figure  6.  One  can  see  Chat,  within  the  limits  of  scaled  distance  covered  by  Dewey 
and  Sperrazza,**  Sachs*  law  is  indeed  verified  by  their  tests. 

In  experiments  for  scaled  distances  close  to  explosive  sources  and  under  very  low  ambient  pressure 
conditions.  Jack  and  Armendt’  showed  that  the  entire  character  of  the  blast  wave  changes  at  altitude 
from  sea  level  conditions,  and  that  Sachs*  scaling  for  pres‘='ures  does  not  apply — the  assumption  of  air 
behaving  as  a perfect  gas  is,  of  course,  invalidated  for  these  tests  close  to  the  blast  source.  An 
anomaly  observed  by  both  Olson,  et  al . and  Jack  and  Armendt^  is  that  this  law  apparently  does  apply 
for  the  rotlectod  impulse  parameter,  even  very  close  to  the  explosive  source.  We  note  here  that  this 
agreement  is  strictly  fortuitous,  and  an  explanation  is  given  later. 

Sachs'  law  is  used  almost  universally  to  predict  effects  of  change  in  the  condition  on  blast 
parameters.  .Most  authors  correctly  identify  the  law  as  due  to  Sachs,  but  some,  such  as  Brode®  and 
Glasstone,’  simply  use  it  with  no  mention  of  its  author. 

Inherent  in  both  liopkinson  and  Sachs'  laws,  in  addition  to  the  assumption  in  Sachs’  law  that  air 
behaves  as  a perfect  gas,  are  the  assumptions  that  gravity  and  viscosity  effects  are  negligible. 

Sachs'  law  includes  Hopkinson  scaling  as  a special  case,  when  tliere  are  no  changes  in  ambient  condi- 
tions between  model  and  prototype  experiments. 

For  worst  case  design,  the  sea  level  relations  in  Figure  2 are  recommended.  In  particular,  the 
Z values  ranging  from  0.1  to  1.0  have  some  experimental  verification. 

KEFLtCTtI)  PKESSURE  COEFFICIENT  VERSUS  ANGLE  OF  INCIDENCE 

Wlien  a detonation  occurs  adjacent  to  and  above  a protective  structure  so  that  no  amplification  of 
the  initial  shock  wave  occurs  between  the  explosion  source  and  the  protective  structure,  then  the  blast 
loadings  acting  on  the  structure  are  free-air-burst  blast  pressures. 

As  the  incident  pressure  wave  moves  radially  away  from  tne  center  of  the  explosion,  it  w*  K come 
in  contact  with  the  structure  and,  upon  contact,  the  initial  wave  pressures  and  impulse  are  r.af'^rced 
and  reflected  (Figure  7).  The  reflected  pressure  pulse  of  Figure  7 is  typical  for  infinite  plane 
reflectors. 


Fig. 7 Prcssurc-limc  variation  for  free  air  burst 
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Wlien  Che  shock  wave  Impinges  on  a surface  oriented  so  Chat  a line  which  describes  the  path  of 
travel  of  the  wave  is  normal  to  the  surface,  then  the  point  of  initial  contact  is  said  to  sustain  the 
maximum  (normal  reflected)  pressure  and  impulse.  The  positive  phase  pressure  and  impulse  patterns  on 
the  structure  vary  with  distance  from  a maximum  at  this  normal  distance  Ra  ^ minimum  (incident 
pressure)  where  the  plane  of  the  structure's  surface  is  perpendicular  to  the  shock  front.  The  peak 
pressures,  impulses,  velocities,  and  other  parameters  of  this  shock  environment  for  a spherical  TNT 
explosion  versus  the  scaled  distance  (Z  = are  given  in  Figure  2. 

The  variation  of  the  pressure  and  impulse  patterns  on  the  surface  between  the  maximum  and  minimum 
values  is  a function  of  the  angle  of  incidence  a.  This  angle  is  formed  by  the  line  which  defines  the 
normal  distance  Ra  between  the  point  of  detonation,  and  the  structure  and  the  line  R which  defines  the 
path  of  shock  propagation  between  the  center  of  the  explosion  and  any  other  point  in  question. 

The  effect  of  the  angle  of  incidence  on  the  peak  reflected  pressures  is  shown  in  Figure  8,  which 
is  a plot  of  the  angle  of  incidence  versus  the  peak  reflected  pressure  coefficient  as  a function  of 
the  peak  incident  pressure.  The  peak  reflected  pressure  is  obtained  by  multiplying  the  peak 

reflected  pressure  coefficient  C^x  the  peak  incident  pressure  Pso-  design  purposes,  the  other 

blast  parameters  except  the  duration  of  the  wave  may  be  taken  as  those  corresponding  to  the  reflected 
pressure  and  are  obtained  from  Figure  2.  The  duration  of  the  blast  wave  corresponds  to  the  duration 
of  the  free  air  pressures. 


ANGLE  OF  INCIDENCE  oc  (DEGREES) 


hg.8  Reflected  pressure  cocfncicm  versus  angle  of  incidence 
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The  air-burst  blast  enviroaraent  is  produced  by  detouattoiis  which  occur  above  the  jiround  surface 
and  at  some  distance  away  from  the  protective  structure  so  that  the  initial  shock  wave,  propagating 
away  from  the  explosion,  impinges  on  the  ground  surface  prior  to  arrival  at  the  structure.  As  the 
blast  wave  continues  to  propagate  outward,  a front  known  as  the  Mach  front  (Figure  9)  is  formed  by  the 
interaction  of  the  initial  wave  (incident  wave)  and  the  reflected  wave,  which  Is  the  result  of  the 
reinforcement  of  the  incident  wave  by  the  ground. 


Fig.d  buisl  blast  environment 


Some  variation  of  the  pressures  over  the  height  of  the  Mach  front  occurs  but,  for  design  purposes, 
this  variation  can  be  neglected  and  the  shock  considered  as  a plane  wave  over  the  full  height  of  the 
front.  The  pressure-time  variation  of  the  Mach  front  is  similar  to  that  of  the  incident  wave  except 
that  the  magnitude  of  the  press'ire  is  somewhat  larger. 

The  height  of  the  Mach  front  Increases  as  the  wave  propagates  away  from  the  center  of  the  detona- 
tion. This  increase  in  height  is  referred  to  as  the  path  of  the  triple  point  and  is  formed  by  the 
intersection  of  the  initial,  reflected,  and  Mach  waves.  A structure  is  subjected  to  a plane  wave 
(uniform  pressure)  when  the  height  of  the  triple  point  exceeds  the  height  of  the  structure. 

If  the  height  of  the  triple  point  does  not  extend  above  the  height  of  the  structure,  then  the 
magnitude  of  the  applied  loads  will  vary  with  the  height  of  the  point  considered.  Above  the  triple 
point,  the  pressure-time  variation  consists  of  an  interaction  of  the  incident  and  reflected  Incident 
wave  pressures  resulting  in  a pressure-time  variation  different  from  that  of  the  Mach  incident  wave 
I pressures.  The  magnitude  of  pressures  above  the  triple  point  Is  smaller  than  that  of  the  Mach  front, 

i In  roost  practical  design  situations,  the  location  of  the  detonation  will  be  far  enough  away  from  the 

; structure  so  as  not  to  produce  this  pressure  variation. 

ESTIMATE  OF  EXTERNAL  BLAST  PRESSURES  ON  STRUCTURES 

The  techniques  that  will  be  described  now  were  developed  for  the  civil  or  structural  engineer. 
Since  the  writer  has  never  seen  these  techniques  applied  Co  aircraft,  the  explanation  will  be  developed 
around  a rectangular  structure.  The  application  to  aircraft  shapes  is  therefore  left  to  the  reader. 

One  advantage  of  this  approach  will  be  that  the  reader,  if  he  uses  these  approaches  in  estimating  blast 
loads  on  aircraft,  will  have  no  false  Illusions  concerning  Che  extent  to  which  these  approximate  rela- 
tions have  been  applied  and  experimentally  verified  against  aircraft  structures. 

The  blast  loading  on  a structure  caused  by  a high-explosive  detonation  is  dependent  upon  several 
factors:  (I)  the  magnitude  of  the  explosion,  (2)  location  of  the  detonation  relative  to  the  structure 

in  question,  and  (3)  the  geometrical  configuration  of  Che  structure. 

The  procedures  presented  here  for  the  determination  of  the  external  blast  loads  on  structures  are 
restricted  to  rectangular  structures. 


FORCKS  ACTINCi  ON  STKUCTURtS 


The  forces  acting  on  a structure  associated  with  a plane  shock  wave  are  dependent  upon  both  the 
peak  values  and  the  pressure-t ime  variation  of  tlie  incident  and  dynamic  pressures  acting  in  the  free 
field.  The  peak  values  of  the  free-field  incident  pressures  for  the  several  blast-loading  categories 
were  treated  previously. 

For  each  pressure  range*  there  is  a particle  or  wind  velocity  associated  with  the  blast  wave  that 
causes  a dynamic  pressure  on  objects  in  the  path  of  the  wave.  In  the  free  field,  these  dynamic  pres- 
sures are  essentially  functions  of  tiie  air  density  and  particle  velocity.  For  typical  conditions, 
standard  relationships  have  been  established  between  the  peak  incident  pressure  and  tlie  peak 
dynamic  pressure  The  magnitude  of  the  peak  dynamic  pressure  is  solely  a function  of  the  peak 

incident  pressure  and,  therefore,  independent  of  the  size  of  the  explosion. 

For  design  purposes,  it  is  necessary  to  establish  the  variation  or  decay  of  both  the  incident  and 
dynamic  pressures  with  time  since  the  effects  on  a structure  subjected  to  a blast  loading  depend  upon 
the  intensity-time  history  of  the  loading  as  well  as  on  the  peak  intensity.  The  idealized  form  of  the 
incident  blast  wave  (Figure  1)  is  characterized  by  an  abrupt  rise  in  pressure  to  a peak  value,  a period 
of  decay  to  ambient  pressure,  and  a period  in  which  the  pressure  drops  below  ambient  (negative  pressure 
phase)  . 

The  rate  of  decay  of  tlie  Incident  and  dynamic  pressures,  after  the  passage  of  the  shock  front,  is 
a function  of  tlie  peak  pressure  and  the  size  of  the  detonation.  For  design  purposes,  the  actual  decay 
of  the  incident  pressure  may  be  approximated  by  the  use  of  an  equivalent  triangular  pressure-time  pulse 
The  actual  duration  of  the  positive  pressure  phase  is  replaced  by  a fictitious  positive  duration  which 
is  expressed  as  a function  of  the  total  positive  Impulse  and  the  peak  pressure: 


The  above  relationship  for  the  equivalent  triangular  pulse  is  applicable  to  the  incident  pressures  as 
well  .18  the  reflected  bl.ist  pressures;  liowever,  in  the  case  of  the  latter,  the  value  of  the  impulse 
used  with  Kq . 1 is  equivalent  to  that  associated  with  the  reflected  wave.  The  fictitious  duration  of 
the  dynamic  pressures  may  be  assumed  to  be  equal  to  that  of  the  incident  pressures. 

For  determining  the  pressure-time  data  for  the  negative  phase,  a similar  procedure  as  used  in  tlie 
evaluation  of  the  positive  phase  nuiy  be  utilized.  The  equivalent  negative  pressure-time  curve  will 
have  a time  of  rise  equal  to  0.125  whereas  the  fictitious  duration  ’.s  given  by  the  triangular 

equivalent  pulse  equ.it  ion. 


where  i and  P are  the  total  impulse  and  peak  pressure  ol  the  negative  phase  and  is  the  actual 

negative  phase  duration. 

ABOVtCiROUNI)  RLCTANGULAR  STRUCTURES 

For  any  given  set  of  free-field  Incident  and  dynamic  pressure  pulses,  tlie  forces  imparted 
to  an  aboveground  structure  can  be  divided  Into  three  general  components:  (1)  the  force  resulting  from 

the  Incident  pressure,  (2)  the  force  resulting  from  the  dynamic  pressure,  and  (3)  tlie  reflected  pressure 
resulting  from  the  shock  Impinging  upo:i  an  interfering  surface.  Tlie  relative  significance  of  each  of 
these  components  is  dependent  upon  the  geometrical  configuration  and  size  of  the  structure,  and  the 
orientation  of  the  structure  relative  to  the  shock  wave. 

The  interaction  of  the  incident  blast  wave  with  an  object  is  a complicated  process.  To  reduce  the 
complex  problem  of  blast  to  reasonable  terms,  it  will  be  assumed  liere  that:  (1)  the  structure  is 

generally  rectangular  in  shape,  (2)  the  incident  pressure  of  interest  is  in  the  order  of  200  psi  or 
less,  and  (3)  the  object  being  loaded  is  in  the  region  of  the  Mach  reflection. 

Front  Surface.  For  a rectangular  structure  at  low  pressure  ranges,  the  variation  of  pressure  with 
time  on  the  side  facing  the  detonation  is  Illustrated  in  Figure  10a.  At  the  noment  the  incident  shock 
front  strikes  the  surface,  the  pressure  is  immediately  raised  from  zerj  to  the  reflected  pressure 
which  is  a function  of  the  incident  pressure  and  the  angle  of  Incidc  j between  the  shock  front  and  the 
structure  face  (Figure  8). 

The  clearing  time  t;;.  required  to  relieve  the  reflected  pressures  Is  represented  as 


3£ 
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(3) 


where  is  the  velocity  ot  the  shock  front  and  S is  equal  to  the  height  of  the  structure  or  one-half 
its  width  Wst  whichever  is  smaller.  The  pressure  P acting  on  the  iront  surface  after  time  ' is  the 
algebraic  sum  of  the  incident  pressure  Pg  drag  pressure  c'/yj. 


Ihe  drag  coefficient  Cq  gives  the  relationship  between  the  dyn«iniic  pressure  and  the  total  t rans  la  t ional 
pressure  In  the  direction  of  the  wind  produced  by  tlie  dynamic  pressure,  and  varies  witli  the  Mach  number 
(or  with  the  Reynold’s  number  at  low  incident  pressures)  and  the  relative  geometry  of  the  structure. 

A value  of  Cq  - I for  the  front  surface  is  considered  adequate  for  the  pressure  ranges  considered  in 
this  manual . 
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Pig. 10  P'font  surface  loading 


At  higher  pressure  r.inges,  the  above  procedure  may  yield  a fictitious  pressure-t  Inc  curve  because 
of  the  extremely  short  pressure  pulse  durations  Involved.  Therefore,  the  pressure-time  curve  construc- 
tion must  be  checked  to  determine  Its  accuracy.  The  comparison  Is  made  by  constructing  a second  curve 
(dotted  triangle  as  Indicated  In  Figure  10a),  using  the  total  reflected  pressure  Impulse  depending 
upon  the  shock  environment.  The  fictitious  duration  °f  the  reflected  wave  Is  calculated  from 

2i 

''r  ■ (5) 


where  is  the  peak  reflected  pressure  previously  used. 
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Wliuljfvfi  t urvt'  ^ivtjs  tlio  smallest  value  ul  the  impulse  (area  under  curve),  use  that  curve  in 
lalvtUatiiig  tlie  wall  loading.  The  reflected  pressure  impulse  includes  the  effects  of  both  the  incident 
and  tlyn.jmic  pressures. 

The  pressure-t irae  curve  assumes  that  the  shock  front  of  the  blast  wave  is  essentially  parallel  to 
the  front  surface  of  the  structure  (plane  wave  shock  front).  In  those  cases  wliere  the  detonation  is 
located  relatively  close  to  the  structure,  the  pressure  and  impulse  patterns  acting  across  the  wall 
surface  will  vary.  In  these  cases  an  average  pressure-t ime  curve  (Figure  iOb)  may  be  constructed  using 
average  values  of  the  loading  acting  over  a length  /)£,*  of  the  structure  wiiere  £)£■  is  equal  to  1 . 3 
is  the  normal  distance  between  the  detonation  and  the  front  surface  of  tiie  structure)  but  not  greater 
than  2i>.  In  this  solution,  both  tl»e  pressures  and  the  duration  are  average  values  of  the  loading 
acting  on  the  wall  between  points  1 and  2. 

As  in  the  case  of  the  plane  wave  frot\t,  the  calculations  for  the  average  loading  on  the  front 
surface  subjected  to  a varying  shock  wave  front  must  consider  the  effect  of  extremely  short  pressure 
pulses.  Here  the  values  of  and  tj-  are  taken  as  average  values  of  the  loading  acting  between  points 
I and  2 . 

Roof  and _S_Uie  Surfaces.  As  the  shock  front  traverses  a structure,  a pressure  is  imparted  to  the 
roof  slab  and  side  surfaces  equal  to  the  incident  pressure  at  a given  time  at  any  specified  point 
reduced  by  a negative  drag  pressure.  The  portion  of  the  surface  loaded  at  a particular  time  is  depend- 
ent upon  the  location  of  tl>e  shock  front  and  the  wave  lengths  (L^  and  L^“)  of  the  positive  and  negative 
pulses . 

To  accurately  determine  the  overall  loading  on  a surface  where  the  span  direction  is  perpendicular 
to  the  shock  front,  a step-by-step  analysis  of  the  wave  propagation  should  be  made.  This  analysis 
would  include  a simultaneous  dynamic  analysis  of  the  stresses  produced  in  the  element  at  any  given  time. 
To  simplify  this  procedure,  an  approximate  method  of  calculating  the  pressure-t ime  history  has  been 
developed  wlwre  tt\e  actual  loading  v^n  the  surface  has  been  replaced  by  an  equivalent  uniform  loading, 
which  will  produce  stresses  in  the  element  similar  to  those  produced  by  the  actual  loading  as  the  blast 
wave  crosses  the  surface.  For  this  analysis  to  be  valid,  it  is  assumed  that  the  reinforcement  on  both 
faces  is  continious  across  the  span  length. 

The  peak  value  of  the  Incident  pressure  decays  and  the  wave  length  increases  as  the  shock  wave 
traverses  the  roof.  As  illustrated  in  Figure  11a,  tiie  maximum  stress  in  the  member  occurs  when  the 
shock  front  is  at  point  d although  the  point  of  maximum  stress  is  located  elsewhere.  The  equivalent 
uniform  loading  versus  time  is  shown  in  Figure  11b  and,  to  simplify  the  calculations , lias  been  taken 
as  a function  of  the  blast  wave  parameters  at  point  b,  the  back  end  of  the  element  being  considered. 

The  equivalent  load  factor  Cy  and  the  blast  wave  location  ratio  D/L  are  obtained  from  Figure  12  as  a 
function  of  the  wave  lengtl»--span  ratio  fhe  pressure  builds  up  linearly  from  the  tine  wlien 

the  blast  wave  reaches  tin  beginning  of  the  element  (point  f)  to  time  when  the  blast  wave  reaches 
point  d.  The  peak  v.»lue  Pq  o!  the  pressure  is  the  sum  of  the  contributions  of  the  equivalent  incident 
and  drag  pressures. 


/’ 
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where  is  the  peak  overpressure  occurring  at  point  b and  corresponds  to  the  value  of 

Following  the  peak,  the  pressure  decays  linearly  to  zero  at  time  + t^^f  where  is  the  time  at  which 
the  blast  wave  reaches  the  end  of  the  element  (point  b)  and  tQf  is  the  fictitious  duration  of  the 
positive  phase.  The  negative  phase  of  the  equivalent  uniform  loading,  if  required,  may  be  taken 
conservatively  as  that  occurring  at  point  b. 


The  drag  coefficient  Co  for  the  roof  and  side  surfaces  is  a function  of  the  peak  dynamic  pressure. 
Recommended  values  are  as  follows: 


Pe^  djrnamic_j;»_^essur_e  Draj^  coet  f ic  ietu 

0-25  psi  -0.40 

25-50  psi  -0.30 

50-130  psi  -0.20 


If  the  span  direction  of  the  member  is  parallel  to  the  siiock  front,  the  pressure-time  curve  shown 
in  Figure  13  should  be  used  where  L is  the  width  of  the  strip  or  element  being  considered. 

Rjjar  Surface.  As  tlie  shock  front  moves  past  the  rear  edges  of  the  roof  and/or  side  surfaces,  the 
pressure  front  will  expand  forming  secondary  waves  which  propagate  over  the  rear  wall.  In  the  case  of 
long  buildings,  this  second.iry  wave  enveloping  tlu*  back  surface  essentially  results  fiom  tne  spillover 
from  the  roof,  wliile  in  short  buildings,  surface  loads  arise  from  the  Interaction  of  the  spillover  both 
from  the  roof  and  the  side  surfaces.  in  both  cases,  the  secondary  waves  are  reinforced  due  to  their 
Impingement  with  reflecting  surfaces.  Presently,  little  information  is  available  on  the  overall  effects 
on  the  rear  surface  loading  produced  by  the  reflections  of  the  secondary  waves. 
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SHOCK  FRONT  LOCATION 
FOR  MAXIMUM  STRESS 


o)  SECTION  THROUGH  STRUCTURE 


b)  AVERAGE  PRESSURE  - TIME  VARIATION 


I'igll  Roof  Jtid  Side  surface  K»ading  (span  duectu»n  perpendicular  U»  shock  tiorri  I 


111  most  design  cases,  the  primary  reason  for  determining  the  blast  loads  acting  on  llie  rear 
surface  is  to  determine  the  overall  drag  effects  (noth  front  and  rear  surtace  loadings)  on  the  struc- 
ture. For  tills  purpose,  a procedure  may  be  used  where  the  blast  loading  on  tlie  surface  is  calculated 
using  the  equivalent  uniform  metiiod  utilized  for  computing  the  blast  loads  on  the  roof  and  side  sur- 
faces. Here  the  peak  pressure  of  tl\e  pressure-time  curve  (Figure  l^b)  is  calcuiaied  using  the  peak 
pressure  tliat  would  occur  at  distance  Hg  (point  o on  Figure  l*4a)  past  the  rear  edge  of  the  roof 
slab.  The  ecjulvalent  load  factor  is  based  on  the  wave  leiigtii  of  tiie  pi-ak  pressure  above  and  tiu* 
unsupported  length  of  the  rear  surface,  as  are  the  time  i*f  rise  and  duration. 

Like  the  roof  and  side  surfaces,  the  blast  loads  acting  on  the  rear  surface  are  a function  ol  the 
drag  pressures  in  addition  to  the  incident  pressures.  The  dynamic  pressure  ot  the  drag  cv'rresponds  to 
that  associated  with  the  equivalent  pressure,  while  the  recommended  drag  ct'e  f 1 ic  i ents  are  tlie  same  as 
used  fi*r  the  roof  and  side  surfaces. 


IMPLLSl  MOMIMtM  SIKUTIRAL  RL.SI*ONSL 

The  classical  approach  to  structural  response  for  pressure  loads  that  are  of  slnirt  time  duration 
with  respect  to  the  lowest  period  of  the  responding  structure  is  called  impulse  momentum  technique, 
riiis  approach  Is  a rational  one  in  that  it  utilizes  Newton's  equations  of  motion. 

Before  discussing  the  fund.imenlal  principles  of  dynamic  analysis,  the  principles  used  in  the 
analysis  of  structxires  under  static  load  will  be  reviewed  briefly.  Two  different  methods  are  used 
either  separately  or  concurrently  In  static  analysis:  one  is  based  on  the  principle  of  equilibrium, 

and  the  other  on  wi>rk  done  and  internal  energy  considerations. 
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{-ig.I2  I qinvuleni  load  tacior  and  blast  v^ave  location  ratio  versus  wave  length-span  ratio 


L’n«lk'r  the  application  of  external  loads,  a given  structure  is  deformed  and  internal  lorces  devel- 
oped in  its  members.  In  order  to  satisfy  static  equi  1 ibt  ium,  tlie  vector  sum  of  all  tlie  external  and 
internal  forces  acting  on  .my  free  body  portion  of  the  structure  must  he  equal  to  ;»ero.  For  the 
equilibrium  of  the  structure  as  a whole,  the  vector  sum  ol  the  external  forces  and  the  reactions  of 
the  foundation  must  also  he  equal  to  zero. 

The  method  based  on  work  done  and  energy  cons i de rat i ons  is  sometimes  used  when  it  is  necessary 
to  determine  tlie  deform.it  ion  of  a structure.  In  this  method,  use  is  made  i>f  the  fact  that  the  deforma- 
t iiin  of  the  structure  c.iuses  tlie  point  of  application  of  the  external  load  to  he  displaced.  The  force 
then  does  work  on  the  stnuture.  Meanwhile,  because  of  the  structural  deformations,  potential  energy 
is  stored  in  the  structure  In  the  fomi  of  strain  eniTgy.  By  the  principle  of  conservation  of  energy, 
the  work  done  by  the  external  force  and  the  energy  stored  in  tlie  roerahers  must  be  equal.  In  static 
analysis,  simplified  methods  such  as  the  method  of  virtual  work.  ni\d  the  method  of  the  unit  lo.ad  arc 
derived  from  the  general  principle  of  energy  conservation. 

In  the  analysis  of  statically  indeterminate  structures,  in  addition  to  satisfying  the  equations 
of  equilibrium.  It  is  necessary  to  include  a calculation  of  the  deformation  of  the  structure  in  order 
to  arrive  at  <i  complete  solution  of  tlie  internal  forces  In  the  structure.  The  methods  based  on  energy 
considerations,  such  as  the  method  of  least  work  and  the  method  based  on  Cas t tgl iano 's  theorems,  are 
generally  used. 


i 
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SHOCK  FRONT 


a)  SECTION  THROUGH  STRUCTURE 


b)  AVERAGE  PRESSURE  - TIME  VARIATION 

Rear  surface  loading 


analysis  of  structures  under  dynamic  loading,  tl.e  same  two  methods  basically  are  used 


(7) 


where 


f - total  external  force  as  a function  of  time 
R = total  Internal  force  as  a function  of  time 
M * total  mass 

«3  » acceleration  of  the  mass. 


As 

kinetic 


for  the  principle  of  conservation  of  energy, 
energy  and  the  strain  energy: 


the  work  done  must  be  equal  to  the  sum  of  the 


y/D  » KE  + SE 


(8) 
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whore 


KE 

SE 


work  done 
kinetic  energy 
strain  energy 


and  the  strain  energy  includes  both  the  reversible  elastic  strain  energy  and  the  irreversible  plastic 
strain  energy.  Thus  the  difference  between  structures  under  static  and  dynamic  loads  is  the  presence 
of  inertial  force  (.Ma)  in  the  equation  of  dynamic  equilibrium,  and  of  kinetic  energy  in  the  equation 
of  energy  conservation.  Both  terras  are  related  to  the  mass  of  the  structure.  Hence,  the  mass  of  the 
structure  becomes  an  important  consideration  in  dynamic  analysis. 

In  the  above  discussion,  the  two  metliods  of  analysis  arc  described  as  if  they  were  independent  of 
each  other.  This  is  not  the  case.  The  force  and  energy  equations  are  interdependent  and  are  derived 
from  each  other.  Although  these  two  methods  are  not  independent,  the  conveniences  of  applying  them  to 
a given  problem  varies  depending  on  tlie  particular  problem  under  consideration.  Generally,  the  force 
equation  is  convenient  for  analysis  of  structures,  and  the  energy  equation  is  convenient  for  design. 

DYNAMICALLY  LQUIVALtNI  SYSTLMS 

In  dynamic  analysis  there  are  only  tliree  quantities  to  be  considered:  (1)  the  work  done,  (2)  the 

strain  energy,  and  (3)  tl)e  kinetic  energy.  To  evaluate  the  work  done,  Llie  displacement  at  any  point 
on  the  structure  under  external  distributed  or  concentrated  loads  must  be  known.  The  strain  energy  is 
equal  to  the  summation  of  the  strain  energies  in  all  the  structural  elements,  wliich  may  be  in  bending, 

compression,  shear,  or  torsion.  The  kinetic  energy  involves  the  energy  of  translation  and  rotation 

of  all  the  masses  of  the  structure.  Tlie  actual  evaluation  of  tliese  quantities  for  a given  structure 
under  dynamic  load  would  be  complicated.  However,  for  practical  problems  this  can  be  avoided  by 
using  appropriate  assumptions. 

In  order  to  simplify  the  problem,  a given  structure  is  replaced  by  a dynamically  equivalent  system. 
The  distributed  masses  of  the  given  structure  are  lumped  together  into  a number  of  concentrated  masses. 
The  strain  energy  is  assumed  to  be  stored  in  several  weiglitless  springs  which  do  not  have  to  behave 
elastically.  Similarly,  the  distributed  load  is  replaced  by  a number  of  concentrated  loads  acting  on 
the  concentrated  masses.  Therefore,  the  equivalent  system  consists  merely  of  a number  of  concentrated 

masses  joined  together  by  weightless  springs  and  subjected  to  concentrated  loads  which  vary  with  time. 

This  concentrated  mass-spring- load  system  is  defined  as  an  equivalent  dynamic  system. 

The  reduction  of  a given  structure  to  an  equivalent  dynamic  system  involves  the  principle  of 
dynamic  similarity  which  is  the  requirement  that  the  work  done,  strain  energy,  and  kinetic  energy  of 
the  equivalent  system  must  be  identical,  respectively,  witli  tliose  of  the  given  structure.  Any  given 
structure,  a multistory  building  or  an  individual  structural  element  such  as  a beam  or  column,  can  be 
approximated  by  an  equivalent  dynamic  system  for  the  purpose  of  dynamic  analysis. 

DYNAMIC  ANALYSIS  AND  DtSICN 


The  analysis  of  a structural  member  under  applied  dynamic  loads  consists  of  the  determination  of 
the  deflection  of  the  member  using  Newton's  equation  of  motion.  Analysis  is  a direct  procedure  while 
design  is  a verification  by  analysis.  The  acceleration,  velocity,  and  displacement  versus  time  rela- 
tionships are  obtainable  from  the  mathematical  solution  of  the  equation  of  motion  of  the  structural 
system.  But  the  solution  also  requires  expression  of  the  load  and  structural  parameters. 

In  the  typical  design  problem,  the  load-time  relationship  is  known  while  the  structural  parameters 
are  not  and,  therefore,  must  be  assumed.  The  displacement  of  the  structural  element  at  certain  critical 
points  can  be  related  to  moments  or  stresses  at  these  same  points.  Solution  of  the  equation  of  motion 
will  give  the  m.ixlmum  deflection  of  the  system  which  can  be  compared  to  acceptable  or  limit  values  and, 
if  adequate,  the  structural  properties  originally  assumed  are  satisfactory  for  the  loading  parameters. 

If  not,  the  assumption  of  a new  set  of  structural  properties  will  give  a new  equation  of  motion  which, 
in  turn,  can  be  solved.  The  process  of  structural  design  is  a t r ial-and-error  procedure — a series  of 
analyses  to  verify  assumed  designs. 

It  can  be  seen  that  the  writing  and  solution  of  the  equation  of  motion  for  a structure  are  essen- 
tial, but  usually  mechanical,  steps  In  the  design  process.  They  are  not  always  readily  apparent  in 
the  equations  and  charts  presented  here,  but  they  are  nevertheless  there. 

The  simplest  dynamic  system  consists  of  a concent ratc'd  mass  supported  by  a weightless  spring  and 
subjected  to  a concentrated  load  as  shown  in  Figure  15.  A single  displacement  variable  .v  Is  sufficient 
to  describe  its  motion;  hence  this  system  Is  called  a slnglo-degree-of-f rcedom  system.  The  degree  of 
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freedom  of  a dynamic  sjystem  is  defined  as  the  number  of  independent 
specify  completely  the  configuration  of  the  system.  This  manual  is 
single-degree~of-freedom  systems  only. 


displacement  variables  needed  to 
restricted  to  consideration  of 


Fig.15  Typjcal  single-degree-of-freedotn  system 


Two  fundamental  methods  are  available  for  treating  simple  systems  subjected  to  dynamic  forces. 

The  first  of  these  methods  is  concerned  with  solving  the  differential  equations  of  the  system  by  either 
classical,  numerical,  or  grapliical  means.  The  second  method  of  analysis,  a chart  solution,  depends 
on  solutions  that  have  been  determined  by  use  of  the  first  method  and  is  an  approximate  solution  to 
the  problem. 

CiLNLRAL 

Structural  elements  must  develop  an  Internal  resistance  sufficient  to  stop  all  motion  within  the 
limits  of  deflection  prescribed  for  the  particular  design;  flexural  design  is  the  determination  of  this 
required  resistance.  The  load  capacity  of  tl»e  member  depends  on  the  peak  strength  developed  by  the 
member  and  on  the  ability  of  the  member  to  sustain  its  resistance  for  a specific  though  relatively 
short  period  of  time.  This  section  describes  tlie  methods  used  for  flexural  design  of  structural 
elements  and  is  divided  Into  two  parts:  (1)  elements  that  respond  to  the  pressure  only,  and  pressure- 

time relationship  that  corresponds  to  the  low-  and  intermediate-pressure  design  ranges,  respectively, 
and  (2)  elements  that  respond  to  the  Impulse  which  fall  Into  the  high-pressure  design  range. 

Both  the  method  of  analysis  and  the  maximum  deflection  criteria  are  different.  Elements  that 

respond  to  the  pressure  only  and  pressure-time  relationship  are  designed  using  a response  chart,  while 

elements  that  respond  to  the  impulse  are  designed  using  either  impulse  charts  for  large  deflection 
designs  or  an  impulse  method  for  designs  with  limited  deflections.  Maximum  allowable  deflections  for 
the  first  case  are  based  on  the  ductility  ratio,  and  in  the  second  case,  on  the  angle  of  rotation  at 
the  supports. 

I LKMhNTS  THAT  Rl  SPONI)  TO  PRtSSURE  ONLY  ANO  PRESSURE-TIME  RELATIONSHIP 

These  elements  were  previously  shown  to  be  those  where  the  true  magnitudes  of  the  pressure  and 
duration  of  the  applied  blast  loads  roust  be  known.  The  determination  of  the  dynamic  response  of  these 

systems  Is  accomplished  using  a response  cliart  which  relates  the  dynamic  properties  of  the  blast  loads 

(pressure  and  duration)  to  those  of  the  element  (natural  period  of  vibration,  resistance,  and  deflection). 


19 


The  maxlmuni  allowable  deflection  is  expressed  in  terms  of  the  ductility  ratio  or  the  ratio  of  the 
maximum  deflection  to  the  equivalent  maximum  elastic  deflection  X£  of  the  element.  Maximum  deflec- 
tions associated  with  the  analysis  of  elements  that  respond  to  pressure  only  and  elements  that  respond 
to  pressure-time  relationship  usually  conform  to  those  designated  as  limited  deflections. 

In  order  to  utilize  these  response  charts,  both  the  blast  load  (pressure-time  history)  and  tue 
resistance-deflection  curve  are  idealized  to  linear  or  bilinear  functions.  The  idealization  of  these 
functions  consists  of  replacing  the  true  functions  by  one  or  a series  of  straight  lines. 

Idealization  of  Load  and  Resistance  Functions.  The  most  consnonly  used  idealized  pressure-t ime 
function  is  an  Initially  peaked  triangle  defined  by  a peak  pressure  B and  a duration  T.  This  simplified 
loading  function  has  been  studied  extensively.  The  basic  consideration  in  replacing  the  given  loading 
by  an  idealized  loading  shape  is  that  the  maximum  displacement  produced  by  each  must  be  equal.  No 
general  rules  can  be  given  for  the  idealization;  h<n*^ever,  the  following  will  serve  as  a guide.  Since 
t;7j  is  the  time  at  which  the  maximum  deflection  Xjn  occurs,  the  idealized  load  and  the  given  load  should 
be  very  close  to  each  other  in  the  time  interval  from  zero  to  t;n  have  the  same  area  under  the 
respective  curves  in  the  same  time  interval. 

For  example,  consider  the  solid  loading  curve  shown  in  Figure  16.  For  a time  to  reach  maximum 
deflection,  the  idealized  loading  curve  is  best  defined  by  B\  and  , while  for  a time  t;n2  ^he  loading 
would  be  defined  by  B2  and  73. 


'I 

TIME 


hg.l6  id(ali/u(ion  ol  pfcssure-limc  functions 


The  typical  resistance  function  for  limited  deflection  design  consists  of  elastic,  elas to-plast ic , 
and  plastic  ranges,  and  for  the  large  majority  of  problems,  can  be  adequately  Idealized  to  a bilinear 
form.  For  the  particular  case  where  the  structural  element  exhibits  only  elastic  and  plastic  behavior, 
idealization  of  the  resistance  function  Is  not  necessary. 

Natural  Period  of  Vibration.  Wlien  designing  for  limited  deflections  using  the  response  chart, 
the  effective  natural  period  of  vibration  is  required.  This  effective  period  of  vibration,  when 
related  to  the  duration  of  a blast  loading  of  given  intensity  and  a given  structural  resistance, 
determines  the  m.ixlmum  transient  deflection  Xf^  of  the  structural  clement. 


The  effective  natural  period  of  vibration  is 


where  m^,  is  the  effective  unit  mass  and  is  the  equivalent  unit  stiffness  of  the  system. 

The  values  used  for  and  K^r  for  a particular  element  depend  upon  the  allowable  maximum  deflec- 
tions permitted.  For  example,  when  designing  for  completely  elastic  behavior,  the  elastic  values  of 
the  effective  unit  mass  and  unit  stiffness  would  be  used;  for  elasto-plastic  action,  the  effective 
unit  mass  is  the  average  of  the  elastic  and  elasto-plastic  values  and  the  stiffness  is  the  equivalent 
unit  stiffness  For  design  in  the  plastic  range,  a weighted  value  of  the  effective  unit  mass  is 

used  with  the  equivalent  unit  stiffness  K^'.  When  plastic  deformations  less  than  tliose  specified  by 
the  limited  deflection  criteria  are  used,  the  use  of  the  average  value  of  the  average  elastic  and 
elasto-plastic  effective  unit  masses  and  the  plastic  effective  unit  mass  will  suffice  for  most  struc- 
tural elements  considered  in  this  manual. 

Determination  of  Maximum  Resj>_o_ns_e . The  maximum  response  of  a pressure  or  pressure-t  ime  sensitive 
system  with  an  elasto-plastic  resistance  function  subjected  to  a triangular  loading  pulse  is  shown  in 
Figure  17.  This  nondimens ional  chart  is  a convenient  plot  of  solutions  made  by  numerical  integration 
of  the  equation  of  riK'tion  using  a digital  computer.  If  any  two  of  the  four  parameters  are  known,  the 
remaining  two  can  be  found.  It  sliould  be  noted  that  since  the  blast  load  is  expressed  in  terms  of 
pressure  (force/unit  area),  unit  values  are  used  for  the  resistance  and  stiffness  of  an  element. 

In  the  typical  design  example,  the  pressure-time  loading  is  calculated,  and  this  curve  is  then 
Idealized  to  the  triangular  form  defined  by  B and  T.  If  a structural  member  is  assumed,  the  idealized 
resistance  function  defined  by  X£-,  and  can  be  determined  along  with  the  natural  period 
Then,  knowing  the  ratios  B/r^  and  T/r.v,  Xm  and  t/n  can  be  obtained  readily  from  Figure  17  by  reading 
the  ratios  Xm/Xfr  and 

Tile  deflection  is  the  maximum  deflection  corresponding  to  the  assumed  parameters.  The  value 
of  is  the  time  at  which  the  maximum  deflection  occurs  and  is  used  to  clieck  the  adequacy  of  the 
idealization  of  the  loading  curve.  If  either  t[„  or  Xu^  are  unsatisfactory,  the  procedure  is  repeated. 

Determination  of  Rebound.  In  the  design  of  elements  that  respond  to  the  pressure  only  and  pressure- 
time relationship,  the  element  must  be  designed  to  resist  the  negative  deflection  or  rebound  which  can 
occur  after  the  maximum  positive  deflection  has  been  reached.  The  ratio  of  the  required  unit  rebound 
resistance  to  the  ultimate  unit  resistance  so  that  tiie  element  will  remain  elastic  during 

rebound,  is  obtained  from  Figure  18.  Entering  with  tlie  ratios  Xjf,/X£  and  T/Ty  previously  determined 
for  the  design,  the  required  unit  rebound  resistance  r~  can  be  read  in  terms  of  the  originally 
designed  ultimate  unit  resistance  r^. 

It  may  be  noted  that  If  the  loading  is  applied  in  a relatively  short  time  compared  to  the  natural 
period  of  vibration,  the  required  rebound  resistance  can  be  equal  to  the  resistance  in  the  initial 
design  direction.  When  the  loading  is  applied  for  a relatively  long  time,  the  maximum  deflection  is 
reached  when  the  positive  forces  are  still  large,  and  the  rebound  resistance  is  reduced. 


fcLtMhNTS  THAT  RtSPONI)  TO  IMPULSE 

General  Equations  for  Maximum  Response.  Wlien  an  element  responds  to  the  impulse,  the  maximum 
response  depends  upon  the  area  under  tlie  pressure-t Ime  curve  (impulse  of  the  blast  loading).  The 
magnitude  and  time  variation  of  the  pressure  are  not  Important.  Response  charts  based  on  pressure- 
time relationships  are  therefore  not  applicable  to  these  problems.  Instead,  the  element  resistance 
required  to  limit  tlie  maximum  deflection  to  a specified  value  is  obtained  through  the  use  of  a semi- 
graphical  method  of  analysis. 

Consider  the  pressure-  and  resistance-time  curves  shtiwn  in  Figure  19.  The  resistance  curve  de- 
picted is  for  a two-way  element  with  a resistance-deflection  function  having  a post-ultimate  range. 

From  .%'ewton's  equation  of  motion  it  can  be  shown  that  the  summation  of  the  areas  (considering  area  A 
as  positive  and  area  B as  negative)  under  the  time  curves  up  to  any  time  divided  by  the  corresponding 
effective  masses  is  equal  to  the  Instantaneous  velocity  at  that  time: 
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The  displacement  at  time  t^  is  found  by  multiplying  each  differential  area  divided  by  the  appropriate 
effective  mass  by  Its  distance  to  t^  and  summing  the  values  algebraically: 
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Fig. 19  Pressure  time  and  resistance-time  curves  for  elements  that  respond  to  impulse 


For  an  element  to  be  In  equilibrium  at  its  maximum  deflection,  its  impulse  capacity  must  be 
numberically  equal  to  the  impulse  of  the  applied  blast  load.  With  the  use  of  the  foregoing  equations, 
the  expressions,  which  define  the  motion  and  capacity  of  elements  subjected  to  impulse  type  loads,  can 
be  defined.  These  expressions  are  presented  for  both  large  and  limited  deflection  criteria. 

Case  1 — Large  Deflections . 

1.  Utilizing  Eq.  11  and  by  taking  moments  of  the  areas  under  the  pressure-  and  resistance-time 
curves  (Figure  19)  about  time  assuming  that  the  unit  blast  impulse  1^,  is  applied  instantaneously 
at  t ® 0,  and  that  the  time  to  reacli  yield  ty  is  also  close  to  zero,  the  expression  for  the  maximum 
deflection  Xm  is 


r t , \ t - -r-l  r 
u 1 \ m 2 / 
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If  moments  of  the  areas  are  taken  about  tj^,  then  the  deflection  at  partial  failure  , is 
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(13) 


Using  Eq.  10  and  summing  the  areas  to  and  recognizing  that  the  instantaneous  velocity  at  equals 
zero 


a 0 


(14) 
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2.  Solution  of  the  above  three  simultaneous  equations  is  accomplished  by  solving  £q . 14  for  t„i 
and  substituting  this  expression  into  Eq.  12,  and  solving  Eq . 13  for  t|  and  substituting  this  expression 
into  the  modified  Eq.  12.  After  combining  and  rearranging  terms,  the  general  response  equation  becomes 
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(15) 


The  left  side  of  this  equation  is  simply  the  Initial  kinetic  energy  resulting  from  the  applied  blast 
impulse  and  the  right  side  is  the  modified  potential  energy  of  the  element.  The  modification  is  re- 
quired since  the  above  analysis  requires  the  use  of  two  equivalent  dynamic  systems  (before  and  after 
time  ti).  The  modification  factor  equates  the  two  dynamic  systems.  If  the  effective  mass  in 

each  range  was  the  same,  would  equal  one  and  the  right  side  of  the  expression  would  be  21rX 

which  is  the  potential  energy. 

3.  For  one-way  elements  which  do  not  exhibit  the  post-ultimate  resistance  range,  or  for  two-way 
panels  where  the  maximum  deflection  Xjj,  is  less  than  L‘q.  13  becomes 


u 


5.  The  above  solutions  are  valid  only  for  large  deflection  design  (support  rotations  >5  degrees) 
since  the  variation  of  resistance  with  deflection  in  the  eiasto-plast ic  range  has  been  ignored.  This 
limitation  Is  based  on  the  assumption  that  the  time  to  reach  yield  ty  and  the  duration  of  the  impulse 
tQ  are  small  in  comparison  to  t;„. 


CRITICAL  IMPULSL  WITHIN  CRITICAL  TIML 

Another  technique  which  has  received  limited  use  is  the  critical  impulse  within  a critical  time. 
An  excerpt  from  two  references  (Reference  10  and  11)  is  herein  Included  which  illustrates  the  method. 
From  Reference  10,  the  following. 


CRniCAL  IMPLLSL 


The  identification  and  selection  ot  the  critical  time  within  which 
an  impulse  must  be  received  by  a target  In  order  to  have  d.unage  inflicted 
is  essentially  empirical.  However,  there  are  available  some  theoretical 
guidelines.  Thus,  for  a simple  systcan  at  rest  but  capable  of  harmonic 
motion,  a maximum  velocity  and  maximum  amplitude  of  swing  are  given  by 
an  Impulse  of  very  short  duration.  Increasing  the  duration  of  the  given 
Impulse  results  In  a decrease  In  achieved  velocity  and  amplitude.  However, 
this  decrease  Is  relatively  small  until  the  duration  of  the  impulse 
becomes  about  one-quarter  of  the  natural  period  of  oscillation,  when 
further  Increases  give  marked  decreases  in  .unplitude.  These  observations 
have  been  amply  confirmed  by  analog  computer  studies  as  illustrated  in 
FIGURE  2.  In  addition,  it  may  be  noted  that  a simple  harmonic  oscillator 
travels  from  zero  to  maximum  displacement  in  one-quarter  of  Its  natural 
period  and  that  the  enforced  motion  of  a structural  member  leading  to 
damage  effects  should  not  be  greatly  different. 

On  a basis  of  these  considerations,  it  appears  that  a time  equal  to 
one-quarter  of  the  natural  period  of  oscillation  should  for  any  struc- 
tural element  represent  a rough  sort  of  cut-off  time  for  impulse  effec- 
tiveness. Representative  calculations  based  on  the  determination  of 
such  natural  periods  are  discussed  in  the  following  section. 


In  addition  to  structural  deformations  as  a form  of  blast  damage, 
there  is  also  the  situation  where  an  object  such  as  a truck  or  a body  is 
physically  displaced  and  hurled  from  its  original  position.  The  forces 
giving  rise  to  such  a hurling  action  must  at  least  be  sufficiently  great 
to  override  the  forces  of  gravity,  and  they  are  arbitrarily  taken  to  be 
those  forces  which  by  acting  for  one  second  can  produce  the  corresponding 
free-fall  velocity.  A time  in  the  order  of  one  second  thus  becomes  the 
critical  time  for  blast  damage  received  by  this  hurling  iDechanism. 
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ratio  impulse  duration 

NATURAL  PERIOD 

FigL’RI-  2 Krlaiivr  amphiudr  of  undamped  harmonic  motion  vnsus  duration  of  applied 
impulse. 


To  sammarize,  the  critical  times  for  explosive  blast  damage  are 
arbitrarily  but  reasonably  taken  as  one-quarter  of  the  natural  period 
for  structures  and  structural  members,  and  one  second  for  destructive 
hurling  away  for  nonattached  objects. 

Critical  Impulse 

Structural  damage  inflicted  by  explosive  blast  is  the  result  of  an 
impulsive  load  that  exceeds  the  resistance  of  a material.  In  ordinary 
impulsive  loading,  a stress  wave  moves  througl);  the  material  such  that 

0 = QCV,  (7) 

where 

0 is  the  stress, 

p is  the  density, 

C is  the  velocity  of  sound  in  the  material,  and 

V is  the  relative  particle  velocity  associated  with  the  stress  wave 
(the  product  pC  is  known  as  the  "acoustic  impedance") 

For  the  limiting  situation  where  the  stresses  developed  are  equal  to 
the  dynamic  yield  strength,  there  also  is  a limiting  or  critical  particle 
velocity,  Vq%  which  if  exceeded  causes  permanent  displacement  within  the 
material . 

Experimental  values  for  this  critical  particle  velocity  are  available 
for  the  tensile  failure  of  some  of  the  common  metals.  For  otiier  materials, 
the  critical  particle  velocity  can  in  principle  be  found  from  the  dynamic 
yield  strength  when  such  dynamic  yield  values  are  known  (they  are 
appreciably  greater  than  the  conventional  static  values).  However, 
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blast  produced  failures  are  seldom  ones  in  simple  tension.  Hence, 
the  critical  particle  velocity  applicable  in  the  disruption  of  a target 
by  blast  is  not  so  simply  described,  but  It  should  be  at  least  closely 
related  to  the  critical  particle  velocity  for  failure  in  simple  tension. 
For  purposes  here,  these  two  critical  particle  velocities  will  be  con- 
sidered as  being  identical. 


For  blast  damage,  the  critical  particle  velocity  is  actually  a dif- 
ferential velocity  referenced  to  a part  of  the  system  at  rest.  Consider, 
for  example,  a structural  system  where  a skin  is  stretched  between  two 
rigid  support  members.  When  an  impulsive  load  is  applied,  the  Inertia 
of  the  support  members  keeps  them  from  moving  significantly,  while  the 
thinner  skin  is  accelerated  to  a velocity  that  may  exceed  the  critical 
velocity,  in  which  case  tearing  of  the  skin  results.  Similar  considera- 
tions apply  in  systems  that  cantilever  out  from  a massive  inertial 
reference  system,  such  as  the  empennage  or  the  wing  of  an  aircraft 
cantilevered  out  from  the  fuselage. 


The  critical  impulse  per  unit  area  associated  with  the  critical 
velocity  is  given  by  elementary  considerat ions  as 


(8) 


where 


p is  the  density  of  the  target  material, 
6 Is  the  thickness  of  the  material,  and 
is  the  critical  velocity. 


The  critical  Impulse  may  also  be  expressed  in  terms  of  the  dynamic 
yield  strength  for  the  material  by  combining  the  above  relations  to  give 


6a  /C 

y 


<9) 


where 


C is  the  velocity  of  sound  in  the  material,  and 


CTy  Is  the  dynamic  yield  strength. 


The  velocity  of  sound,  C,  can  be  replaced  by  its  equivalent 


/E/P 


where 


£ is  the  modulus  of  elasticity 

and  this  gives  a third  form  for  the  expression  for  critical  Impulse 


(10) 


These  three  forms  of  the  expression  for  critical  Impulse  are  all 
expressed  In  coherent  units  and  are  all  equivalent.  The  first  of  these 
would  appeal  to  those  who  have  worked  with  impulsive  loads  and  shock 
waves  In  metals,  but  the  others  may  on  occasion  be  preferred  since  they 
contain  items  directly  related  to  the  ordinary  strength  characteristics 
of  a material . 


The  above  expressions  apply  as  part  of  the  criterion  for  blast 
damage  when  the  damage  occurs  by  deformations  of  a structure  or  a 
structural  member.  For  the  Important  situation  where  a target  Is 
damaged  by  hurling  It  bodily  from  Its  original  position,  the  critical 
impulse  Is  arbitrarily  taken  to  be  that  which  could  cause  the  entire 
object  to  reach  free-fall  velocity  in  one  second.  This  is  convenient 
mathematically,  and  also  Is  In  «iccord  with  the  observation  that  the 
corresponding  fall  of  about  16  feet  is  ordinarily  very  damaging  to  most 
targets  of  Interest. 


RtPRtStNTATIVt  CALCULATIONS 


Values  for  the  blast  damage  criterion,  critical  Impulse  within  a 
critical  time,  were  computed  both  for  a parked  aircraft  and  for  a light 
Industrial  building.  For  the  aircraft,  three  types  of  damage  were 
considered:  (1)  tearing  the  skin,  (2)  breaking  off  a wing  section,  and 

(3)  physical  displacement  by  hurling  the  entire  plane. 


Aircraft  Skin  Failure 


For  damage  by  the  mechanism  of  tearing  the  skin,  the  critical  time 
Is  taken  from  a computed  value  for  the  resonant  frequency  of  this  struc- 
tural member.  The  skin  was  assumed  to  be  of  1/16-lnch-thlck  aluminum, 
clamped  to  stringers  on  8-inch  separation.  For  clamped  panels  of 
infinite  length,  the  resonant  frequency,  /,  is  given  by  the  relation 

f = 217,600  X 0.985  x 6/l>^  (11) 

where  the  factor  217,600  pertains  to  this  type  of  structure,  and  the 
factor  0.985  Is  for  the  type  of  metal. 

Here 

6 is  the  skin  thickness  in  inches,  and 
h the  distance  of  spacing  In  Inches. 

Thus , 

f = 217,600  X 0.985  x (1/16)(8)^  = 212  cps  (12) 

This  resonant  frequency  corresponds  to  a natural  period  of  vibration  of 
1,000/212  = 4.70  milliseconds  (ms),  and  to  a critical  time  of  4.70/4  = 1.2  ms 

The  critical  impulse  required  for  tearing  the  skin  is  computed  in 
any  one  of  three  separate  ways,  depending  on  the  primary  data  available. 

Data  on  the  aluminum  skin,  as  obtained  from  various  sources  are  as  follows. 

5 (thickness)  * 1/16  inch  *=  0.0052  feet 

P (density)  = 2.7  g/cm^  = 2.7  x 62.4/32.2 

= 5.22  slugs/ft  3 

C (speed  of  sound)  * 16,470  fps 

Vq  (critical  velocity)  = 240  fps 

O-y  (dynamic  yield  strength)  = 140,000  psl  » 2.02  x 10^  psf 

E (elastic  modulus)  “ 11  x 10^  psi  = 1.58  x 10^  psf 

1.  Critical  impulse  from  critical  velocity: 

I = p6v  * 5.22  X 0.0052  x 240  = 6.5  psf-sec 
c c 

= 45  psi-ms  (13) 

2.  Critical  impulse  from  dynamic  yield  strength  and  speed  of  sound: 

= 6py/C  = 0.0052  X 2.02  x 10^/16,470  - 6.4  psf-sec 
* 44  psi-ms  (14) 

3.  Critical  Impulse  from  elastic  modulus  and  dynamic  yield  strength: 

I * (p/£)*^^6nv  = (5.22/1.58  X 10“*)^^^  x 0.0052  x 2.02  x 10^ 
c 

= 6.0  psf-sec  = 42  psi-ms  (15) 

Differences  in  the  calculated  results,  which  in  this  instance  are  minor, 
are  caused  by  disparities  in  the  original  data. 

A diagram  Illustrating  skin  rupture  as  caused  by  a blast  impulse 
that  exceeds  the  critical  Impulse  for  the  skin  within  the  critical  time 
is  shown  schematically  in  FIGURE  3. 

The  computed  minimum  overpressure  for  this  damage  becomes  43/1.2 
= 36  psi,  which  must  persist  for  at  least  as  long  as  1.2  ms  in  order  to 
provide  the  necessary  minimum  impulse. 


Fiouxf.  3.  Rupture  of  Mruciural  skin  panel.  C'^mputed  critical  impulse.  43  psi-ms;  com- 
puted critical  time.  1.2  ms. 
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An  excerpt  from  Reference  11  gives  background  to  support  the  ability  to  estimate  critical  fracture 
stress  and  critical  impact  velocity. 


Scabbing — When  the  primary  longitudinal  wave,  which  is  one  of  com- 
pression, strikes  the  lower  edge  of  the  plate  it  is  reflected  and  comes 
back  as  a wave  of  tension.  If  the  wave  is  decaying  in  intensity, 
interference  between  the  Incident  and  reflected  waves  will  cause  a tensile 
stress  to  be  built  up  a short  distance  from  the  bottom  of  the  plate. 

The  tension  may  reach  a high  enougli  value  to  cause  the  plate  to  fracture. 
Wlien  a fracture  of  this  type  occurs,  the  plate  is  said  to  have  scabbed. 
Fracturing  of  this  type  was  first  described  by  Hopkinson  many  years  ago 
and  is  sometimes  known  as  a hopkinson  fracture. 

The  mechanics  of  scabbing  can  be  treated  qualitatively  in  a fairly 
elementary  fashion.  Consider  a laterally  infinite  slab  of  finite  thickness 


(a)  (b) 


F'iir.  9-.1  Stress  Distribution  Near  Free  Sur- 
face of  riate  of  Finite  Thicknt'ss  and  Infinite 
l.aleral  F.xtent.  Uefore  refleetum  ami  O'V 

shortly  after  r«'fleeti«m  «»f  transient  disturt»ance. 


and  assume  a perfectly  elastic  plate  material.  Suppose  that  the  pressure 
pulse  Is  of  the  form  shown  in  Fig.  9-M;  that  Is,  the  pressure  rises 
Instantaneovisly  and  then  begins  to  decay.  Wlien  the  disturbance  strikes 
the  free  (right)  boundary  of  the  plate,  it  will  be  reflected  as  a tension 
wave.  The  incident  compression  wave  and  its  reflected  tension  counter- 
part will  interfere.  The  resultant  distribution  of  stress  within  the 
plate  at  a particular  time  during  reflection  of  the  wave  will  be  that 
shown  in  Fig.  9- 3b.  The  tension  AB  will  increase  as  the  reflected  wave 
moves  to  the  left.  If  at  some  point  the  metal  can  no  longer  support 
the  tension,  it  will  fracture  and  a scab  will  fly  off.  For  a body  with 
a definite  critical  normal  fracture  stress,  the  thickness  6 of  the  scab 
will  be  equal  to  one-half  the  distance  within  the  wave  that  corresponds 
to  a decrease  In  stress  equal  to  the  critical  normal  fracture  stress. 

Studies  of  the  scabbing  of  metal  plates  have  shown  that  the  two 
factors  that  influence  scabbing  most  are,  (a)  the  shape  of  the  stress 
wave,  and  (b)  a critical  normal  fracture  stress  that  Is  characteristic 
of  the  material  acted  upon.  The  latter  may  be  influenced  by  the  con- 
ditions of  loading  and  the  state  of  stress  in  the  body. 

The  very  strong  dependence  that  the  location  of  the  scab  has  on 
the  shape  of  the  stress  wave  is  perhaps  best  illustrated  by  some  examples 
of  scabbing  of  flat,  mild  steel  plates  that  were  explosively  loaded  In 
the  manner  shown  in  Fig.  5- 3c  of  Chapter  5.  The  plates  in  the  test 
ranged  In  thickness  from  1/4  to  2 inch;  the  explosive  charges  were  1/8 
and  1/4  Inch  in  thlcknes.s.  Results  are  shown  In  Fig.  9-4.  The  most 
interesting  things  shown  by  the  two  curves  are  that  the  thickness  of 
Stab  increases  rather  than  decreases  with  plate  thickness  and  that, 
under  st»me  conditions,  a 1/8  Inch  layer  of  explosive  can  produce  a 
thicker  scab  than  a 1/4  inch  layer.  The  general  qualitative  aspects 
of  the  evirves  of  Fig.  9-4  are  in  agreement  with  expected  behavior. 


l.e.,  flattening  and  elongation,  of  the  pulse  as  it  moves  through  the 
plate.  Decrease  in  steepness  of  the  stress-time  curve  behind  the  wave 
front  with  Increasing  plate  thickness  will  cause  increasingly  thicker 
scabs  to  form.  The  pulse  will  be  more  sustained  and  the  initial  pressure 
somewhat  higher  for  the  thicker  charge;  hence,  the  thickness  of  the  scab 
generated  by  the  1/4  inch  layer  of  explosive  will  be  greater  than  that 
generated  by  the  1/8  inch  layer  when  the  plate  is  thin. 


Plate  Thickness,  inches 


Fi|t.  9-4  -Dependence  of  Scab 
Thickness  on  Plate  Thickness  for  Two 
Thicknesses  of  Explosive.  Plates  were 
mild  steel  and  loaded  in  manner  of 
Fiif.  5-3c. 

Experimentally  determined  values  of  the  critical  normal  fracture 
stresses  for  five  metals  are  listed  in  Table  9-1.  The  relationship 

O = p cv  Equation  (9-1) 

can  be  used  to  associate  a critical  differential  particle  velocity 
with  each  critical  normal  fracture  stress.  Values  of  critical  particle 
velocities  obtained  in  this  way  are  also  listed  in  the  same  table. 

There  is  some,  but  not  conclusive,  evidence  that  indicates  that  the 
state  of  stress  that  exists  in  the  region  of  the  fracture  at  the  time 
the  fracture  occurs  may  affect  the  value  of 


Table  9-1 


Critical  Normal  Fracture  Stresses  and  Associated 
Critical  Impact  Velocities 


Critical  normal 

Associated 

fracture  stress 

critical  impact 

velocity 

psi 

( ft/sec) 

24S-T4  Aluminum 

140,000 

202 

Copper 

410,000 

264 

Brass 

310,000 

216 

1020  Steel 

160,000 

84 

4 1 30  Steel 

440,000 

2 35 

AIR  BLAST  INSIDKCOMPARTMKNTS 


Warhead  Lffcc! 

The  subject  of  air  blast  inside  the  various  compartments  of  the  aircraft  such  as  the  wing  box, 
fuselage  fuel  tanks,  and  the  various  fuselage  compartments  is  usually  restricted  to  the  smaller  warheads. 
A rule  of  thumb  sometimes  used  is  that  a 2-pound-weight  charge  placed  inside  an  aircraft  will  destroy 
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the  aircraft.  This  rule,  ol  course,  has  been  found  experimentally  to  have  exceptions  for  certain  air- 
craft compartments  where  the  local  geometry  by  shielding,  angle  of  Incidence  effects,  and  venting  can 
increase  the  amount  of  explosive  required  to  kill  the  aircraft  by  almost  a factor  of  10. 

b'ffects  lo  Consider 

1.  The  warhead  length-to-diameter  ratio  will  have  a very  pronounced  effect  for  the  close-scaled 
distances  in  many  aircraft  compartments.  It  can  be  highly  directional. 

2.  Momentum  as  well  as  impulse  must  be  considered. 

3.  The  fragmentation  can  be  quite  different  within  five  diameters  of  the  warhead.  (This  can 
produce  a fragment  hole  pattern  which  the  wounded  structure  must  respond  to  in  its  maximum  structural 
deflect  ion . ) 

Nothing  is  available  to  predict  the  effects  of  a contact  burst  (where  the  warhead  is  in  con- 
tact, or  partially  through  the  structure,  at  detonation). 

5.  Detonator  location,  multiple  detonators,  possible  jacket  focusing  effect,  etc.,  can  have  a 
pronounced  influence  on  effective  explosive  yield  in  various  angular  directions. 

b.  Data  on  the  ability  of  aircraft  fasteners  to  resist  tension  created  by  blast  loads  are 
nonexistent.  Hence,  predictions  of  amount  of  aircraft  skin  removal  (due  to  internal  blast)  is  difficult 
if  not  impossible. 

Best  Method  of  PrtK'edure 

Despite  all  the  uncertainties  previously  listed,  we  must  proceed  and  do  tl»e  best  we  can  if  we  are 
ever  to  arrive  at  a rational  approach.  The  following  approach  is  tentatively  suggested: 

1.  Consider  reflection  from  adjacent  walls  (using  weight  of  explosive  (WK)  + 202). 

2.  Check  to  see  whether  volume  is  vented. 

3.  Proceed  with  structural  response  calculation. 


Keflectii>n.s  From  Adjacent  Walls  of  Closed  Vented  Compartments 

In  confined  areas,  amplification  of  the  Initial  shock  front  (free-air  pressures)  due  to  reflections 
within  the  structure  will  occur.  At  any  given  point  on  a particular  reflecting  surtace,  the  total 
impulse  loading  is  a combination  of  the  contributions  from  the  initial  shock  and  from  tlie  shock  re- 
flected from  adjacent  surfaces. 

An  approximate  method  for  the  calculation  of  the  total  impulse  In  a fully  vented  cubicle  has  been 
developed  in  connection  with  the  Safety  Design  Criteria  Program.  The  total  reflected  Impulse  acting 
on  various  points  of  each  surface  of  the  cubicle  have  been  calculated  and  then  Integrated  to  obtain  the 
total  Impulse  loads.  In  order  to  simplify  the  calculation  of  the  response  of  a cubicle  wall  to  these 
applied  loads,  the  total  Impulse  Is  assumed  to  be  distributed  uniformly,  giving  an  average  value  of  the 
impulse  acting  on  any  one  surface.  The  actual  distribution  of  the  blast  impulses  is  highly  irregular 
due  to  the  multiple  reflections  and  time  phasing,  and  results  in  localized  high  shear  stresses  in  the 
elements.  When  designing,  the  use  of  average  Impulse  loads  Is  predicted  on  the  ability  of  tlie  element 
to  transfer  these  localized  loads  to  regions  of  lower  stress. 

The  method  of  calculating  the  average  blast  impulses  was  developed  using  a theoretical  procedure 
based  on  semiempi rlcal  blast  data  and  on  the  results  of  response  tests  of  reinforced  concrete  slabs. 

The  calculated  average  impulse  loads  have  been  compared  with  those  obtained  from  the  results  of  tests 
of  a scale-model  steel  cubicle  and  have  shown  good  agreement  for  a wide  range  of  cubicle  arrangements. 

The  parameters  which  are  necessary  to  determine  the  average  impulse  loads  are  structure  configura- 
tion and  size,  charge  weight,  and  charge  location.  Surfaces  depicted  are  not  frangible  for  determining 
the  shock  loadings  but  could  be  frangible  when  considering  gas  pressure  buildup.  For  the  four-wall 
cubicle  with  roof,  one  of  the  surfaces  shown  must  be  frangible  or  vented  for  such  gas  pressures. 

Figure  20  shows  the  definition  of  a four-wall  cubicle. 

Because  of  the  wide  range  of  required  parameters,  the  procedure  for  the  determination  of  the 
impulse  loads  was  programmed  for  solution  on  a digital  computer.  The  results  of  these  calculations, 
presented  in  Figures  21  through  29,  give  the  scaled  average  unit  blast  impulse  Ifc,  as  a function  of  the 
parameters  defining  the  configurations  presented  in  Figure  20.  Each  illustration  is  for  a particular 
combination  of  values  of  h/H , 1/L,  and  the  number  N of  reflecting  surfaces  adjacent  to  the  surface  for 
which  the  Impulses  are  being  calculated.  The  wall  (if  any)  parallel  to  the  surface  in  question  was 
found  to  have  a negligible  contribution  to  the  Impulse  for  the  range  of  parameters  used  and  was  therefore 
not  considered. 
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NOTES: 


1 B DENOTES  BACK  WALL.  S SIDE  WALL  AND  B ftOOf 

2 NUMBERS  <N  PARENTHESES  INDICATE  NUMBER  N OF  REFLECTING 
SURFACES  ADJACENT  TO  SURFACE  IN  QUESTION 

i h IS  ALWAYS  MEASURED  TO  THE  NEAREST  REFLECTING  SURFACE 

4 £lS  ALWAYS  MEASURED  TO  THE  NEAREST  REFLECTING  SURFACE 

EXCEPT  FOR  THE  CANTILEVER  WALL  WHERE  IT  IS  MEASURED 
TO  THE  NEAREST  FREE  EDGE 


i FOR  VALUES  OF  SCALED  AVERAGE  UNIT  BLAST  IMPULSE 

BARRIER  AND  CUBICLE  ARRANGEMENTS  SHOWN  SEE  FIGURES  21 
THROUGH  FIGURE  73  REQUIRED  PARAMETERS  ARE 
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Fig. 20  Four-wall  cubicle  with  roof 

The  general  procedure  for  use  of  the  above  Illustrations  is  as  follows: 

1.  Determine  the  values  of  the  parameters  Indicated  for  the  selected  surface  of  the  barrier  or 
cubicle  configuration  and  calculate  the  following  quantities: 


h.  L L L. 

H * L * W * * 

A 


and 


A 


.1/3 


2.  In  most  cases  the  above  procedure  will  require  interpolation  for  one  or  more  of  the  parameters 
which  define  a given  situation,  in  order  to  obtain  the  correct  average  impulse  load. 

3.  Because  of  the  limitation  in  the  range  of  the  test  data  and  the  limited  number  of  values  of 
the  parameters  given  in  the  above  impulse  charts,  extrapolation  of  the  data  given  in  Figures  21  throu^. 

29  may  be  required  for  some  of  the  parameters  Involved.  On  the  other  hand,  the  limiting  values  as  given 
in  the  charts  for  other  parameters  will  not  r€*quire  extrapolation.  The  values  of  the  average  impulse 
loads  corresponding  to  the  values  of  the  parameters  which  exceed  their  limiting  values  (as  defined  by 
the  charts),  will  be  approximately  equal  to  those  corresponding  to  the  limiting  values.  The  following 
are  rec«»mraended  procedures  that  will  be  applicable  in  most  cases  for  either  extrapolation,  or  establish- 
ing the  limits  of  impulse  loads  corresponding  to  values  of  the  various  parameters  that  exceed  the 
limits  of  the  charts: 

(a)  Plot  curve  of  values  of  versus  for  constant  values  of  L/R^,  L/H , h/H , and  1/L. 
Extrapolate  curve  to  Include  value  of  corresponding  to  the  value  of  2^  required. 

(b)  Extrapolate  given  curve  for  constant  values  of  L/H,  h/H,  and  I/L,  to  include  v.ilue  of 
corresponding  to  the  value  of  L/R^  required. 

(c)  Values  of  corresponding  to  values  of  L/H  greater  than  b shall  be  taken  as  equal  to 

those  corresponding  to  L/H  • 6 for  actual  values  of  h/H,  and  I/L,  but  with  a fictitious  value  of 
L/R/^  in  which  R^  Is  the  actual  value  and  L is  fictitious  value  equal  to 

(d)  Values  of  corresponding  to  values  of  1/L  less  than  0.10  and  greater  than  0.75  shall  be 

taken  as  equal  to  those  corresponding  to  1/L  * 0.10  and  0.75,  respectively. 

(e)  Values  of  corresponding  to  values  of  h/H  less  than  0.15  and  greater  than  0.75  shall 

be  taken  as  equal  to  those  corresponding  to  h/H  » 0.15  and  0.75,  respectively. 
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Fig. 27  Scaled  average  unit  blast  imp'iKe  (N  = 4, 


I mn  mu  ms  u ti  n n Bsttittut  mt;  ou: 
inainB8mtiitnnnsattttisu:mc 
iBomanttnuBitmisutmessi 

Sn8BBttiiaBifnimtttt»n8:, 

mm  liimiHyiffiimrHil 


mnewnniiiimnHB 

mHiHHgiliiiHiiHg 


SStkSM 


iiSi£itt!iiii£HBi;miiing 


imaiiaiRagBMainBmi.r'nnRKWk'aiBniB: 

iwBiiuBaaiWBaiBawnBnawniiiuiBaaa 

UnaHPiBRRBBBmBnHBflRIItnBRBn 

lUiiiiiiiHHmiiiilltiinHilttilMHIl 

ragiiimBiiyimi 

iiiiiiiiinwwHiiiiiiiiiiiiyiiiitiiBiiiii 

liiUiimiHiiliiRllliiimilHHgii 

imiiniiiiygiiimii 


lO^OO- 


F ig.’S  Scaled  average  utnt  hlast  impulse  (N  - 4.  I/L  = 0.-5  and  0 75.  h H 0 ^0l 


40 


L^LianaBi 


mmm 


iii3nWk»ixn^^^ 

inBii.'TUgSiifaraaiiMgB.taniBaimBMil 
iiiiiiiiii  j^BjgK^^Ksarnai'aMMBBW 


ig.2‘>  Scaled  avejagc  unit  blast  impulse  IN  4 


41 


As  previously  mentioned*  a protective  element  subjected  to  the  high-pressure  range  effects  of  an 
explosion  may  be  designed  for  the  impulse  rather  than  for  the  pressure  pulse  only  if  the  duration  of 
the  applied  pressures  acting  on  the  entire  element  is  short  in  comparison  to  its  response  time.  ‘An 
estimate  of  this  duration  may  be  obtained  by  adding  the  time  increments  correspoi.  .ing  to  the  time 
required  for  the  blast  wave  to  fully  engulf  the  element  and  the  duration  of  the  blast  load  at  the 
section  of  the  element  furthest  renK)ved  from  the  explosion.  This  relationship  is  represented  by: 


t 

o 


(‘A  ^ 


(17) 


where 


to  » duration  of  load  (ms) 

• time  of  arrival  of  the  blast  wave  at  the  point  on  the  element  furthest  from  the 
explosion;  defined  by  the  Largest  slant  distance. 

• time  of  arrival  of  the  blast  wave  at  the  point  on  the  element  nearest  to  the 
explosion;  defined  by  the  normal  distance. 

(to);-  * duration  of  the  blast  pressure  at  the  point  on  the  element  furthest  from  the  explosion. 

The  time  of  arrival  of  the  blast  wave  for  the  two  points  of  interest  as  well  as  the  duration  of  the  load 
at  the  furthest  point  on  the  element  are  obtained  from  Figure  2.  To  account  for  the  increased  load 
duration  due  to  the  multiple  reflections  of  the  blast  wave  within  the  cubicle*  the  value  of  has 

been  increased  by  50%. 

A protective  element  subjected  to  the  high-pressure  range  effects  of  an  explosion  may  have  to  be 
designed  for  the  pressure  pulse  rather  than  for  the  impulse  effects  if  the  duration  of  the  applied 
pressures  acting  on  the  entire  element  Is  long  in  comparison  to  the  response  time  of  tlie  element 
(tp,  less  chan  3to) . For  these  cases*  the  actual  pressure-time  relationship  resulting  from  a pressure 
distribution  on  the  element,  which  is  lughly  irregular  due  to  the  multiple  reflections  and  time  phasing, 
may  be  approximated  by  a fictitious  peaked  triangular  pressure  pulse.  The  duration  of  this  pulse 
is  obtained  from  Eq.  17  and  a fictitious  peak  pressure  Pf  is  established  as  a function  of  this  duration 
and  the  average  blast  impulse  i^,  (Figures  21  through  29)  acting  on  the  element. 


The  above  solution  for  the  average  impulse  load  does  not  account  for  increased  blast  effects 
produced  by  contact  charges. 


Partially  Vented  bxplosions 

Interior  Blast  Loadings.  When  the  openings  in  a cubicle-type  structure  are  small  compared  to  the 
total  surface  area*  full  venting  of  a contained  explosion  will  not  occur.  The  effect  of  the  reduction 
of  the  venting  areas  from  a fully  vented  configuration  is  an  increase  In  the  interior  pressures, 

Impulse*  and  duration.  In  addition  to  the  initial  shock  wave  and  its  reflections,  the  accumulation 
of  gases  from  the  explosion  results  in  a pressure  buildup*  and  since  the  venting  area  is  limited,  the 
duration  of  the  pressure  is  prolonged. 

A typical  pressure-time  record  at  a point  on  the  surface  of  a partially  vented  chamber  is  shown 
in  Figure  30.  The  high  peaks  are  the  multiple  reflections  of  the  initial  shock  wave.  The  mean  pressure, 
denoted  as  P/y,*  has  a long  duration  when  compared  to  the  duration  of  the  shock  wave  pressures.  The  maxi- 
mum mean  pressure  is  used  as  the  basis  for  design  and  is  a function  of  the  charge  weight*  contained 

volume  of  the  chamber*  and  venting  area. 

Figure  31  shows  an  experimentally  fitted  curve  based  upon  test  results  of  partially  vented  chambers 
with  small  venting  areas.  The  maximum  mean  pressure  is  plotted  against  the  charge-volume  ratio.  The 
sizes  of  the  venting  areas  used  In  the  individual  tests  were  relatively  small  compared  to  the  volume 
of  the  chambers  and  the  quantities  of  explosive*  and  had  no  noticeable  Influence  on  the  maximum  pres- 
sures recorded. 

Since  little  data  are  available  concerning  the  duration  of  these  pressures*  it  is  suggested  that* 
for  purposes  of  design*  the  pressures  be  assumed  to  have  a very  long  duration  and  a constant  magnitude 
insofar  as  the  response  of  the  structure  is  concerned.  The  total  impulse  used  in  the  design  of  such  a 
cubicle  Is  a combination  of  the  contributions  from  the  initial  shock  wave,  its  reflections*  and  the 
pressures  caused  by  the  gas  accumulations.  Because  of  the  long  duration  aspect  of  the  gas  accumulation 
pressures,  the  strength  of  the  individual  elements  of  the  structure  must  be  equal  to,  or  greater  than* 
the  pressure  Pmo*  whereas  the  response  of  the  elements  to  the  Initial  blast  wave  and  its  reflections 
will  be  similar  to  the  response  of  cubicles  where  full  venting  Is  achieved. 


4: 


Leakage  Pressure  Loadings.  Little  data  are  currently  available  for  the  de tenni nat ion  of  leakage 
pressures  from  a partially  vented  explosion.  Table  2 (from  Reference  12)  lists  the  experimental  venting 
or  "leak"  areas.  In  most  cases,  rather  small  leakage  areas  will  vent  the  structure. 


Kig.30  Pressure-nme  varuMun  for  a pariiall>  venleU  explosion 
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Pig.3l  Test  results  from  four  chambers,  pressure  versus  charge-volume  ratio 


TABLI  2 

Ihmcnsions  of  the  Chambers  t'sed  in  the  Tests  and  Their  Leak  Areas 


Mark 

Form 

Dimensions 

Volume,  m^ 

Leak  area, 

A 

Sphere 

Dlam  ■ 0.92  m 

O.Al 

15-80 

B 

Tube 

D 1 am  » 1.5  ra 
Length  - 2.5  m 

4.5 

30-400 

C 

"Cube" 

2.5  X 2.3  X 2.3  m 

13.2 

25-200 

D 

Tube 

Diara  * 2.0  m 
Length  ■ 8.5  m 

28.2  ' 
1 

7-2000 

L 
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SUMMARY  AND  LIMlTAnONS 

1.  The  Impulse,  peak  pressure,  and  durations  of  blast  waves  from  spherical  charges  in  air  are 
listed. 


(a)  Cylindrical  charge  deviations  from  spherical  behavior  are  covered. 

(b)  Momentum  and  traveling-charge  effects  are  mentioned,  although  no  quantitative  effects 
are  listed. 


(c)  Hopklnson  versus  Sachs'  spherical  charge  scaling  laws  (to  account  for  altitude  variations) 
are  discussed. 

2.  The  effect  of  angle  of  blast  wave  Incidence  on  lowering  of  reflected  Impulse  is  discussed. 

3.  Simple  methods  of  computing  external  blast  loads  on  rectangular  structures  are  discussed. 

■4.  A sunanary  of  the  single-degree-of-f reedom  structural  response  is  listed.  Suggestions  as  to 
how  to  handle  the  elastic-plastic  response  range  are  included.  j 

5.  Critical  Impulse  within  a critical  time  is  dicussed.  \ 


I 


b.  Techniques  (assuming  rigid  surface  reflections)  used  to  estimate  internal  blast  loads  are 
covered. 


7. 

cussed. 


Pressure  buildup  due  to  products  of  combustion  in  partially  vented  closed  chambers  is  dls- 
Limited  experimental  data  are  shtn.m. 
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HYDRAULIC  RAM  PRESSURE  PREDICTION 
AND  STRUCTURAL  RESPONSE 


INTRODUCTION 

During  penetration  of  an  aircraft  fuel  cell.  Inillets  and  otlier  liigli-speed  projectiles 
generate  intense  pressure  waves  in  the  liiel.  Response  ol  the  fuel  cell  walls  to  these 
pressure  waves  can  be  catastrophic  failure  due  to  severe  fracturing  of  entrance  and  exit 
fuel  cell  walls.  This  failure  is  often  accentuated,  in  brittle  aluniinuni  cell  walls.  b\  the 
bullet  or  high-speed  projectile  wound  in  the  fuel  cell  wall,  and  a catastrophic  brittle 
fracture  failure  at  low  stress  levels,  can  occur.  This  ithenoinenon.  termed  ’‘hydraulic 
ram.”  is  of  particular  importance  to  the  siir\ ivabilits  of  mihtarc  aircraft. 

A computer  program  was  developed  b\  the  Naval  Weapons  Center  which 
calculates  the  pressure  waves  generated  b>  decelerating,  tumbling  projectiles  in  a fluid. 
The  theory  of  hydraulic  ram  is  given  in  Reference  1.  1 he  pressure  prediction  theorv  in 
this  report  is  excerpted  from  Reference  I.  I he  ludr.iuhc  ram  structural  response 
theory  is  described  later  in  this  report. 

THEORETICAL  CONSIDERATIONS 


To  calculate  fluid  pressures,  it  is  necessary  to  know  the  velocity,  rate  of  kinetic  energy 
loss,  and  time  of  projectile  arrival  as  functions  of  distance  along  the  trajectory.  Estimates  of 
these  quantities  can  be  obtained  by  using  a simple  model  of  the  bullet  behavior. 

The  bullet  is  assumed  to  travel  in  a straight  line,  and  its  deceleration  is  described  by 
Newton’s  Second  Law. 


m 


(1) 


where 

m = bullet  mass 
V = bullet  velocity 
t = time 
D = drag  force 

The  drag  force  can  be  expressed  as 


D = 2 PV-A^,C|3 


(2) 


* Naval  Weapons  Cenler.  hluid  Dvmmu  Aualvus  oj  Hydraulic  Ram  !\\  hy  I A.  l.iindstrom  and  \\  K. 
l ung.  China  Lake.  CA.  NWC.  October  1976.  ( JTC(J/AS-74  I -01  S.  publication  liNCL.ASSH  ILn.l 
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where 

p = fluid  ilensity 
V = bullet  velocity 
Cj)  = drag  coefficient  of  the  bullet 
Aq  = projected  frontal  area  of  an  unyawed  bullet 


Noting  that 


V 


dxb 

dt 


(3) 


where  = bullet  position  along  the  trajectory,  equations  1 and  2 can  be  combined,  yielding 


dV 

dx5 


= -0Vb 


(4) 


Where  the  velocity  decay  coefficient,  /?,  is  defined  as 

the  rate  of  kinetic  energy  loss,  dh/dxb,  where  E = l/2mV3  can  be  expressed  as 

dE  ,,  dV 

dxb  dxb 


(5) 


(6) 


Combining  equations  4 and  6 yields 

^=m/3V2  (7) 

dxb 

Eor  tumbling  bullets,/?  is  a function  of  Xb- 

In  this  model  the  bullet  is  presumed  to  enter  the  test  cell  with  0 degree  yaw  and 
continue  in  this  attitude  with  a constant  drag  coefficient  until  it  reaches  a distance,  X|, 
where  it  begins  to  tumble.  The  bullet  becomes  fully  tumbled  at  a distance,  xt, and  continues 
in  this  attitude  with  a constant  drag  coefficient.  However,  as  evidenced  from  experimental 
results,^  the  bullet  will  impact  the  cell  with  0 degree  yaw  and  continue  to  tumble  along  its 
trajectory  for  a number  of  cycles  before  assuming  a stable  attitude.  The  model  has  an  option 
to  include  this  continuous  tumbling  of  the  bullet.  The  coefficients,  /?j  and  associated 

with  the  0-degree  yaw  and  tumbleil  attitudes,  and  the  value  /?3  is  associated  with  the 
stern-first  attitude.  For  simplicity  it  was  assumed  that  the  tumbling  proceeds  at  a constant 
rate  along  the  trajectory  (Figure  1 );  tliat  is 

^^2  ■ ’‘I  “ ’‘3  ■ ^<2  “ ^4  ‘ X3  = • • • (8) 


^Naval  Weapons  Center.  I'liiiJ  Dynamic  Analysis  of  Hydraulic  Ram  III  by  1'..  A.  Lundstrom  and  W.  K.  l ung. 
China  Lake.  CA.  NWC,  Oelolwr  1974.  (J  IC(:/AS-74 -t-OI  5.  publualion  UNCLASSII  li;i).) 
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Figure  1 . Bullet  Orientation  Versus  Distance  Along 
the  Bullet  Trajectory. 


While  the  bullet  is  tumbling,  i.e.,  xj,  <Xb,  <^2>  varies  radically  as  a function  of  xjj. 
This  variation  can  be  expressed  in  general  form  as 


0 = /3l  +(/32-/Jl)f(Y) 


(9) 


where 


Y = 


Xb-X] 

X2-X1 


(10) 


while 


f(0)  = 0 
f(l)=  1 


An  empirical  function  of  the  following  form  was  used. 

f(Y)=  [l/2- l/2cos(Y)]"  (11) 

The  exponent,  n = 3,  is  varied  (see  Footnote  2)  yielding  the  following  equations,  for 
XI  <Xb<X2 


1/2-  1/2  cos 


Xb-Xl 

X2-X1 


(12a) 


for  X2  < xjj  < X3 
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P(^b^  = P2  (P3  -P2^ 


1/2-  1/2  cos 


Xb-X2 

X3-X2 


3 


(12b) 


When  P is  constant,  equations  3 and  4 can  be  integrated  directly.  In  the  region  0 < x^ 
< xj,  where  P = P\,  the  initial  conditions  are  V = Vq  and  t = 0.  Tlie  integration  then  yields 


Vb  = VoC 


•l^ixb 


and 


‘b-p- 


J_  J_ 

Vb'  Vo 


and  equation  7 yields 


dE 

dxb 


= m/3,Vb- 


(13) 


(14) 


(15) 


In  the  region  of  tumbling,  X]  < Xb  < xt,  algebraic  expressions  cannot  be  found; 
therefore,  a numerical  integration  method  is  used  for  the  model. 

For  stripping  of  bullet  jackets  from  the  AP  (armor  piercing)  core  of  the  API  (armor- 
piercing incendiary)  ammunition,  a crude  method  for  incorporating  the  jacket  energy 
deposition  was  developed.  The  projectile  penetrates  the  fluid  for  a distance,  Xj,  where  the 
jacket  strips.  The  kinetic  energy  of  the  jacket  and  incendiary  material  are  calculated  at  that 
point.  The  energy  deposition  of  the  AP  core  is  calculated  in  the  normal  manner  except  that 
values  of  P appropriate  to  the  core  must  be  used.  The  energy  deposition  of  the  jacket  is 
assumed  to  be  exponential  and  is  added  to  that  of  the  core.  The  equation  for  total  energy 
deposition  is 


dxb 


= niAVb"^  + C 
n 


(16) 


Ejs  = 1/2  mj  Vs2 


(17) 


where 

mj  = mass  of  jacket 

Vj  = velocity  of  bullet  at  stripping  location 
m^  = mass  of  core 

a = a constant  found  to  be  1/3  (Footnote  2) 

Because  of  the  mathematical  difficulties  introduced  by  the  wall  boundary  conditions,  a 
simple  model  (Figure  2)  is  used  which  neglects  the  wall  effects. 
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Figure  2.  Model  of  Drag  Phase  of  Hydraulic  Ram. 


The  bullet  shown  in  Figure  2 is  initially  stationary  in  an  infinite  body  of  fluid  until 
t = 0.  Then  the  bullet  is  impulsively  accelerated  to  an  initial  velocity  of  Vq.  At  times,  t > 0, 
the  bullet  moves  with  a velocity,  V,  in  a straiglit  line  along  the  axis.  It  is  assumed  that  the 
flow  field  can  be  described  in  terms  of  a potential  function,  0,  which  satisfies  the  wave 
equation. 


where  c = speed  of  sound  in  the  fluid. 

Then,  the  fluid  velocity,  u,  is  expressed  as 

iT  = 

and  the  Bernoulli  equation  yields 
P = Po-p  ^-l/2pu2 


at2 


(18) 


(19) 


(20) 


where 


P = total  fluid  pressure 
Pq  = ambient  pressure 
p = fluid  density 
u = I u 

The  boundary  conditions  for  equation  18  are  that  the  fluid  velocity  is  tangential  to  the 
projectile  surface  and  that  P = P,-  on  the  cavity  surface,  where  Pj-  denotes  the  pressure  in  the 
cavity.  It  is  assumed  that  P^  is  constant  throughout  the  cavity;  hence,  the  problem  is  to 
determine  the  pressure  as  a function  of  time  at  any  arbitrary  point  (x,co)  where,  as  shown  in 
Figure  2,  co  is  the  perpendicular  distance  of  this  point  from  the  x axis. 


The  problem  is  further  simplified  by  ignoring  the  boundary  conditions  and  then 
approximating  the  effect  o*"  the  bullet  and  cavity  on  the  fluid  by  the  action  of  a line  of 
sources  distributed  along  the  bullet  trajectory.  Then,  the  resulting  flow  field  is  symmetric 
about  the  x axis,  and  the  potential  due  to  these  sources  can  be  expressed  as 


1 


o 


where 

{ = distance  along  the  trajectory 
r = distance  between  points  £ and  (x,co) 
f = source  strength  at  f 


For  the  finite  sound  speed,  the  integral  must  be  evaluated  along  the  line  t + r/c  = 
constant,  where  t is  the  retarded  time  given  by  r = t - r/c.  Tlie  integration  path  is  included  in 
the  time  space  representation  shown  in  Figure  3.  Using  equation  14,  the  time  of  bullet 
arrival,  tb,  on  the  trajectory  as  a function  of  the  bullet  position,  Xb,  can  be  determined. 
However,  the  results  of  the  theory  are  not  considered  to  be  valid  during  the  cavity  collapse; 
therefore,  the  lower  limit  in  this  integral  can  be  taken  as  zero.  The  upper  limit,  Xb,  denotes 
the  projectile  position  when  t = tb- 


The  source  strength,  f,  is  estimated  by  a method  based  on  the  conservation  of  energy. 
It  is  assumed  that  the  flow  is  confined  to  a section,  dx,  as  shown  in  Figure  4. 


Figure  3.  Integration  Path  of  Integral  in  Equation  21 
on  the  Time-Space  Plane. 


Figure  4.  Flow  Model  for  Estimating  Cavity  Growth  and  Source 
Strength  Variation. 


The  fluid  velocity  ih  tliis  section  is  then 


u = 2iM> 

CO 

At  the  cavity  radius,  co  = a and  u = da/dt  so 


(22) 

(23) 


Following  Birkoffs  theory,'^  the  kinetic  energy,  dk,  of  the  lluid  in  this  section  within  a 
radius,  12,  is 


dk  = 


trp 


'y 

ir^codco 


dx 


(24) 


with  equation  22  and  integrating 

dK  = 4;rpNf-dx  (25) 

where 

N = en(12/a) 

Since  the  upper  limit,  12,  cannot  be  infinite  due  to  the  physical  impossibility  of  allowing  dk 
to  be  infinite,  12  will  be  assumed  to  be  finite.  Because  the  value  of  N varies  slowly  as  a 
function  of  12,  it  can  be  treated  as  a constant.  This  is  justified  for  the  special  case  in  which  a 
bullet  is  traveling  with  a constant  velocity,  since  the  correct  cavity  shape  is  obtained  for 
constant  values  of  the  quantity  12/a  in  the  15  to  30  range.  Physically,  this  step  can  be 
rationalized  by  considering  the  neglected  influence  of  the  noncylindrical  divergence  of  the 
flow. 

The  work,  dW,  done  by  the  difference  between  the  ambient  and  cavity  pressure  is 

dW  = jr(Po  - P(.)a2dx  (26) 

The  energy,  de,  deposited  by  the  projectile  in  the  fluid  at  dx  is 

(t  HirkoO  .iikI  I n /-.ir.inloiu'tlo.  .A /.v.  h'cAc.v.  jnJ  Caiims.  New  York,  Ptcns.  1957. 


with  equations  15  and  16  showing  the  expression  for  dli/dxjj  as  a function  of  the  distance 
along  the  trajectory. 


The  conservatioii  of  energy  can  be  expressed  by  the  relation 


de  = dK  + dW  (28) 


and  with  equations  25,  26,  and  27,  it  becomes 

dx  = 47rpf-Ndx  + TrfP^j  - P^.)a“dx 

(29) 

Defining 

i ^ ,r(Po  - I’e) 

(30) 

® “ pN 

(31) 

and  simplifying,  the  energy  balance  yields  the  expression 

f = ±l/2B\/A2-a2 

(32) 

Eliminating  f from  equation  32,  with  equation  23 

da  .n  \l  a2  2 
a:Tr=  iB  V A“  - a" 
dt 

(33) 

With  boundary  condition  a = 0,  at  the  time  of  projectile  arrival,  tj,, 
33  yields 

integration  of  equation 

± VA2-a2  = A-B(t-tb) 

(34) 

where,  according  to  equation  32 

l/2[BA-B2(t-t,,)] 

(35) 

The  cavity  behavior  reaches  a maximum  radius  when  da/dt  = 0.  From  equation  33,  this 
radius  is  A,  and  from  equation  34,  the  maximum  radius  occurs  at  time,  tp^,  at 

*m  ~ *b  ■* 

(36) 

Tile  effects  of  the  walls  were  neglected,  so  equations  30  and  31  give  only  the  upper  boumls 
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for  A ;itul  B;  licncc  the  source  strentitli,  etiiuilion  35,  will  be  increasintily  inaccurate  as  the 
cavity  approaches  its  maximum  inside  an  actual  luel  cell. 

Now,  the  pressure  field  resulting  from  the  line  of  sources  can  be  calculated.  Substi- 
tuting eciuation  35  for  the  source  strength  in  equation  21  yields 

Xb(r) 

0(x,cj,t)  = -l/2B  J 1 a(?)-B  t-I-tb«)  I df  (37) 

o J) 

Note  that  tlie  retarded  time  t in  equation  21  was  replaced  by  t - r/c. 

To  compute  the  pressure,  the  terms  in  Bernoulli's  equation  (equation  20)  must  be 
evaluated.  Tlie  term,  90/3t,  in  Bernoulli’s  equation  can  be  expressed  in  a simple  form  by 
using  Leibnitz’s  rule  for  the  differentiation  of  integrals.  Equation  37  becomes  (see 
Appendix  A) 

. 1 ^b(r) 


" J o 


where  Rj,  is  the  distance  between  the  bullet  and  the  point  (x,cj)  and  Aj,  denotes  the  value 
of  A,  evaluated  at  Xj,.  The  chain  rule  for  differentiation  gives 

3t  3t  9t  ^ ^ 

where 


where  V is  the  bullet  velocity  evaluated  at  X5,  and  since 
3r_.  1 ^ 1 x-Xt,  9X|,(r) 

3t  c dt  c Rjj  3t 

equation  39  becomes 


1 , 1 9X„(r) 


Solving  for  3X[^(T)/3t  gives 


aXh(r) 

3t 


V X - Xb 
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With  c<iuation  43  and  solving  the  integral,  equation  38  beeomes 


d<t> 


BAu 


^ = .\j2 — 

^ D 


at 


R^. 


\-Xy 


+ l/2B2gn 


X Rq 

X - Xjj  + 


(44) 


Rk 


where  Rq  is  the  distance  between  the  impact  point  and  the  point  (x,co).  In  the  same  manner 
tlie  fluid  velocity  components,  in  the  x direction  and  in  the  co  direction,  can  be 
derived  as 


V X - Xj, 


+ 1/2B 


BAh 

c Rb 

, 

V X - Xb 

■ c Rb 

Xb(r) 

r 

I 

o 

A(£)-B[t- 

BAb 

X ±L 

c Rb 

Rb 

V X - X| 

■ c Rb 

> 

1 j 

' J 

o 

A(£)-B[t- 

d^ 


(45) 


(46) 


It  should  be  noted  that  functions  A($)  and  t|^,(f)  depend  on  the  tumbling  behavior  of  the 
bullet.  Hence,  the  integrals  in  equations  45  and  46  cannot  be  evaluated  explicitly.  Noting 
the  reaction  that 


= u,.^  + u ^ 

A OJ 


(47) 


the  pressure  field  can  be  obtained  by  substituting  equations  44,  45,  46,  and  47  into 
Bernoulli’s  equation  20. 


Waves  reflected  from  rigid  plane  walls  can  be  calculated  exactly  by  means  of  the 
method  of  images.  The  boundary  condition  for  a rigid  wall  is  that  the  normal  component  of 
the  fluid  velocity,  Up,  vanishes  at  the  wall.  Pressure  generated  by  the  bullet  has  been 
approximated  by  a line  of  sources  (or  sinks).  The  boundary  condition  at  the  wall  can  be 
satisfied  by  adding  the  pressure  due  to  a mirror  image  line  of  sources  (or  sinks)  as  illustrated 
in  Figure  5. 


image 


The  boundary  condition  for  a free  surface  is  P = 0 at  the  surface.  A negative  mirror 
c satisfies  this  condition  to  the  extent  that 


J 


55 


» 


1/2u,2 


(48) 


where  uj  and  Uf,  are  the  tangential  and  normal  fluid  velocity  components  at  the  surface. 
Equation  48  is  satisfied  for  most  conditions  during  the  initial  phase  of  hydraulic  ram. 


The  method  of  images  is  easily  extended  to  calculate  reflections  from  the  walls  of 
rectangular  volumes.  The  result  is  a three-dimensional  rectangular  array  of  images.  The 
number  of  images  in  the  array  is  determined  by  the  number  of  wave  refiections  to  be 
included  in  the  calculation.  A method  for  automatic  generation  of  the  image  array  co- 
ordinates was  developed,  and  its  use  is  discussed  in  Airpendix  H. 

The  walls  of  aircraft  fuel  cells  are  neither  free  nor  rigid.  However,  due  to  their  typically 
light  construction,  waves  refiected  from  these  walls  can  be  approximated  by  means  of  the 
methoil  of  images  for  free  surfaces,  fhe  effect  of  the  inertial  properties  of  the  fuel  cell  walls 
on  wave  refiections  is  discussed  in  the  Lundstrom  report  (l-'ootnote  1 ). 


Figure  5.  Geometry  of  Method  of  Images  for  a Single 
Rigid  Surface. 


SUMMARY 

The  pressure  is  calculated  according  to  Bernoulli’s  equation 
P-Po  = -p|f-l/2pu2 

Tlie  first  term  is  dominant  far  from  the  bullet,  while  the  second  term  is  important  closer  to 
the  bullet.  Tlie  boundary  condition  P = 0 at  the  surface  is  satisfied  by  the  method  of  images 
for  the  p 30/3t  term  only.  It  is  suggested  that  the  full  equation  be  used  to  calculate  incident 
pressure  waves  on  the  wall. 

Although  the  action  of  the  cavity  behind  the  bullet  is  accounted  for  in  the  fluid  model, 
the  absence  of  fluid  within  the  cavity  is  ignored.  Therefore,  in  the  absence  of  external 
surfaces,  the  calculated  pressure  will  be  positive  outside  the  cavity,  zero  on  the  cavity 
surface,  and  negative  within  the  cavity.  Indeed,  the  pressure  will  go  to  minus  infinity  as  the 
cavity  axis  is  approached.  In  addition,  the  cavity  presents  an  additional  surface  for  retlecting 
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pressure  w;ives.  Tluis,  waves  arising  from  tlie  bullet,  rellecting  from  tlie  fluid  volume  surface, 
and  then  rcjlcciing  from  ilw  caiily  surface  are  not  accounted  for. 

The  presence  of  the  bounding  free  surfaces  of  the  fluid  volume  often  produces  large 
negative  pressures  within  the  fluid,  thus  producing  bulk  cavitation.  These  phenomena  are 
not  included  in  the  fluid  model. 

It  is  suggested  that  the  pressure  calculated  by  the  program  be  truncated  at  negative 
pressures;  that  is,  if 

P:^0 

then 

P = 0 

This  process  assumed  that  bulk  cavitation  always  occurred  and  also  automatically  accounted 
for  the  absence  of  tluid  in  the  cavity. 

Since  the  bullet  is  modeled  by  a line  of  sources,  the  pressure  goes  to  infinity  at  the 
bullet.  It  is  suggested  that  the  pressure  calculated  by  the  program  bo  truncated  at  the 
stagnation  pressure  of  the  bullet. 


HYDRAULIC  RAM  STRUCTURAL  RESPONSE 


INTRODUCTION 

The  purpose  of  this  section  of  this  report  is  to  describe  modifi- 
cations to  the  methodology  of  the  Hydraulic  Ram  Pressure  Program 
discussed  in  the  previous  section.  These  modifications  allow  for  the 
determination  of  the  wall  deflections  that  result  from  projectile  inter- 
action with  a fluid-filled  rectangular  tank.  The  theory  and  example 
discussed  are  for  the  exit  tank  wall.  However,  the  theory  should  be 
easily  extended  to  the  other  tank  walls  as  required. 

In  the  development  of  the  theory  for  Hydraulic  Ram  Pressure 
Prediction,  the  time  variation  of  the  pressure  can  be  calculated  at 
various  points  internal  to  a fuel  tank  that  has  been  impacted  by  a 
projectile.  The  pressure  prediction  theory  stipulated  that  wall 
reflections  could  be  modeled  by  line  sources  which  were  negative 
mirror  images  (sinks)  of  the  sources  modeling  the  original  projectile 
penetration.  Under  these  conditions,  the  sides  of  the  tank  are 
represented  as  free  surfaces,  and  hence,  by  definition,  the  tank  wall 
pressure  is  zero  at  all  times.  The  hydraulic  ram  structural  response 
theory,  now  to  be  described,  modifies  this  so  as  to  allow  calculation 
of  wall  response  by  the  so-called  "variable  image"  theory. 

STRUCTURAL  RESPONSE  THEORY 


The  specific  response  model  that  is  used  replaces  the  mirror  image 
from  the  exit  wall  used  in  version  one  of  the  code  with  an  infinite  void. 
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Consequently,  a pressure  distribution  p (xyt)  and  a fluid  velocity 
^^(xyt)  can  be  calculated  at  the  position  of  the  exit  wall.  The  struc- 
tural response  of  the  exit  wall  which  acts  against  this  pressure  is  based 
on  Foppls'  calculation.**  The  assumption  is  that  the  total  response  pres- 
sure, q,  generated  by  the  wall  is  such  that 

q = 

where  q-j  is  due  to  the  load  resisted  by  bending  (theory  of  small 
deflections)  and  q^  is  due  to  the  load  resisted  by  membrane  stress.  For 
a given  displacement,  z^,  one  can  calculate  the  total  strain  energy,  U, 
which  is  given  by  the  expression 

U = Uj  + U2 


where 

Uj  = bending  strain  energy 

U2  = membrane  strain  energy 

Timoshenko  et  al**’  ^have  determined  the  analytical  forms  for  the 
membrane  and  bending  strain  energies  under  uniform  pressure  for  a plate 
of  horizontal  length,  a,  and  vertical  length,  b.  The  wall  deflection  is 
assumed  to  have  the  form 


where  x,  y are  the  transverse  coordinates  on  the  wall;  a,  b are  the  asso- 
ciated widths;  and  z^  is  the  center  deflection. 

The  resulting  energy  expressions  are  given  by 


^ 8(12)(1  - J)  \a^  bV 


(49) 


S Itmoshcriko  .tiul  S \k> -Kricgcr  rhtutrv  oj  Plates  and  Shells.  2nd  ltd.  New  York,  McGraw-Hill 

Hook  Compatu.  Iiu  . pp  415  420 

S.  I un»>Nhcnko.  ct  a)  I ihration  f*rohlems  in  f npiu  cntig.  4th  I d.  New  York.  Wiley  and  Company.  1974. 
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where 

V = Poisson's  ratio 
E = modulus  of  elasticity 
h = plate  thickness 
and 


r 2 2 

C = TT 

0 


^ f)  ^ ^ I 

|4(32)(1  + v)  + 9(9  - ^ j 


(bl) 


for 


^2  = 0 
3C  ^ 


The  kinetic  energy  term  for  the  wall  deflection  is  given  by 


T = ill  ab  -2 

' 2g  4 ^0 


where 


vh 

= mass/unit  area  of  plate 


This  allows  the  establishment  of  the  Lagrangian  L = T - U and  its 
associated  differential  equation  so  that 


(at)  (azj  - (szj 


Q = 0 


where  Q represents  the  pressure  forcing  function  due  to  the  projectile- 
fluid  interaction  and  is  related  to  the  total  weighted  average  pressure 
at  the  wall  such  that 


' (!)  "‘’Pt 


where 

b/2  . a/2 

'-b/2  •'-a/2 

From  the  Lagrangian  formulation, 


2 r b/2  ^ a/2 

T = (i)  ^ (?)  (?) 

k / 9 O 


(52) 
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+ ^ (coe)^ 


H 


where 


coe  = TT 


(6-18v)(l.l).24(^.^) 
4(32)  (1  + v)  + 9(9  - v)  TT^  (I  + 


(53) 


(53a) 


The  pressure  distribution  is  calculated  from  the  pressures  due  to  the 
fluid-bullet  interaction  evaluated  at  the  exit  surface.  A problem  with 
this  formulation  is  that  values  are  generated  for  the  wall  velocity 
which  exceed  the  particle  velocity  of  the  fluid.  To  correct  this 
deficiency,  a variable  image  technique  is  employed.  This  allows  the 
interaction  to  be  modified  through  the  fluid  potential  function  in  such 
a manner  as  to  preserve  the  continuity  between  the  plate  velocity  and 
particle  velocity  of  the  fluid. 

Using  the  potential  function,  0^-,  associated  with  the  incident  pres- 
sure wave  and  Bernoulli's  equation,  then 


-p 


!ii  1 2 

at  ■ 2 '"^l 


(54) 


where  ({).  includes  the  potential  produced  by  the  projectile  and  by  the 

waves  reflected  from  the  walls  of  the  container  exclusive  of  the  exit 
wall,  p^  is  the  incident  pressure,  p is  the  density  of  the  fluid  and 

u^  is  the  fluid  velocity  where  u^.  = . Now  assume  that  the  potential 

of  the  reflected  wave  from  the  exit  wall  has  the  form  * = k*.  where  k 

varies  in  time  and  space  and  is  the  mirror  image  of  4)^.. 


Bernoulli's  equation  for  the  reflected  wave  gives 


P 


r 


-p 


(55) 


and  with  similar  substitution 
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dd> 


(56) 


2 

where  n indicates  the  normal  at  the  exit  surface.  Since  the  u term 
in  equation  2 can  be  considered  small  relative  to  the  potential  term, 
it  will  be  ignored  in  these  and  subsequent  equations.  At  the  wall, 
equals  41^  and  34i^/3n  equals  -a‘$'^./an.  From  piston  theory,  a plane 

reflected  wave  is  assumed,  propagating  parallel  to  the  wall  normal  such 
that 


Pr  " "^^‘^rn  ^i  " ^‘^‘^in 

where  c is  the  velocity  of  sound  in  the  fluid. 

Substituting  equation  56  into  equation  57, 

Pr  ■ - ""♦i  (Is)  ■ "P1  ■ (Is)  <“> 

but  from  equation  55, 


Comparing  equations  58  and  59, 


Thus  the  total  pressure  at  the  wall  is 


p^  = p.  + p^  = (1  + k)p.  - (1^) 


(61) 


Similarly,  if  equation  56  is  substituted  into  an  expression  that  equates 
the  wall  and  fluid  velocities  at  the  wall, 


z = u. 


in 


+ u 


rn 


(1  - k)u„, 


62 


or 


Eliminating  tne  derivative  from  p-j., 

Py  = (1  + k)p.  + pcj^d'  - k)u^.  - zj 


To  derive  an  expression  equivalent  to  the  weighted  average  pressure  at 
the  exit  wall,  this  equation  must  be  integrated  by  substituting  into  the 
expression  for  the  forcing  function.  This  gives 


Py  = (1  + k)p.  + pc(l  - k)u^.  - pc^j)  Zq  (62) 


and,  equivalently. 


where  k is  assumed  uniform  at  the  exit  wall.  Given  the  initial 
pressures  and  fluid  velocities  at  specific  time  increments,  one  can 
then  calculate  the  pressure  p.j.  and  the  time  variation  of  k,  z , and  z^. 

These  values  are  then  used  to  predict  a new  set  of  variables  at  the 
next  time  increment  which  are  finally  averaged  with  the  initial  set  to 
yield  the  final  parameters  for  the  given  time  interval.  The  pressures 
used  in  the  calculations  are  gauge  pressures  and  are  allowed  to  go 
negative.  For  time  intervals  exceeding  the  time  for  maximum  deflection, 
this  formulation  occasionally  generates  unrealistic  results.  However, 
experimental  data  and  data  obtained  from  these  calculations  tend  to  be 
reasonably  consistent  for  the  time  intervals  of  interest. 

The  magnitude  of  the  stress  on  the  exit  plate  due  to  these  deflec- 
tions can  also  be  estimated.  From  Hooks  Law,  an  empirical  relation  is 
established  between  the  components  of  stress  and  strain  in  a given 
body.  For  a thin  plate,®  the  equations  relating  the  unit  elongations 

G and  G and  the  normal  stresses  o and  o become 
X y X y 


^ S.  Timo^hc'nko  and  J.  N (toodicr.  Theory  of  lAastieity.  New  York.  McGraw-Hill  Hook  Company.  Inc.. 

1951.  r 
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e 


X 


(64) 


where  E is  the  modulus  of  elasticity  and  v is  Poisson's  ratio.  Using 
the  equations  for  the  plate  displacements  and  the  basic  relations  for 
the  unit  elongations  (e  and  e ) given  by  Timoshenko,'*  at  the  center 

X j 

of  the  exit  plate. 


Ex  = (coe)  ^ and  = (coe)  ^ (65) 


where  coe  is  determined  in  equation  53a.  Solving  equation  64  for 
and  and  substi tuting  equation  65, 


These  are  the  expressions  for  the  stress  used  in  the  code. 

All  these  relations  apply  to  the  special  case  of  a planar  wave  nor- 
mally incident  on  a planar  wall  and  are  assumed  to  be  adequate  for  all 
situations  input  to  the  code.  The  more  general  expression  for  the  pres- 
sure introduced  in  this  section  reduces  to  the  piston  theory  expres- 
sions for  the  pressure  for  transmissive,  rigid,  and  free  surfaces  when 
the  parameter,  k,  has  the  values  of  0,  1,  and  -1,  respectively.  The 
parameter,  k,  is  listed  in  the  output  and  its  variation  models  these 
various  wall  response  conditions. 


EXPERIMENTS 

Experiments  were  performed  with  a 5-ft-high  by  5-ft-wide  by  2-ft- 
thick  water-filled  test  cell.  A number  of  12.7-mm  API  rounds  were 
fired  through  the  2-ft  thickness.  The  exit  wall  was  a 20-inch-square, 
0.125-inch-thick  aluminum  sheet  that  was  rigidly  clamped  at  the  edges. 
Deflection  of  the  exit  wall  was  measured  by  two  high-speed  camera  viewing 
the  exit  wall  profile  from  the  top  and  from  the  side.  Details  of  the 
experiment  are  described  in  a paper  by  Ankeney.^ 


' I).  P.  Ankcncy.  "Hydraulic  Kam  Structural  Response.**  presented  at  the  llvdrodynainic  Ram  Seminar. 
Dayton,  Ohio.  20  22  October  147b. 
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An  equation  of  motion  for  the  exit  wall  deflection  is  previously 
derived  in  this  report.  The  derivation  assumes  a wall  deformation  in 
the  form 


Z = Zq  cos  cos  Y-  (66) 

where  x and  y are  transverse  coordinates  on  the  wall,  a is  the  width  of 
the  square  panel,  and  Zq  is  the  panel  center  deflection. 

Using  the  Lagrange  equation,  a second-order,  nonlinear,  ordinary 
differential  equation  for  Zq  was  obtained.  The  effects  of  plate  mass, 
and  elastic  bending  and  membrane  stresses  are  included.  The  forcing 
function  for  this  equation  is  the  weighted  average  pressure  at  the  wall 


Pt 


(67) 


Assuming  k is  uniform  across  the  plate,  repeating  equation  62  gives 


p.|.  = (1  + k)  p.  + p c(l  - k)  u^.  - p c(J)^  Z 


(68) 


Similarly,  an  equation  for  k is  obtained  by  integrating  equation  63 


9k  _ c 
3t  ■ T. 


2 . 


(1  - - p c(J)  Z 


(69) 


The  weighted  average  incident  pressure  p-j  and  normal  fluid  velocity 
Uni  are  obtained  from  the  previously  discussed  theory.  In  practice,  the 
integrals  for  p-j  and  uni  were  evaluated  by  calculating  p-j  and  Up-j  at 
400  equally  spaced  points  on  the  exit  wall.  Within  the  bullet  cavity 
pi  = -13.8  psi  and  uin  = 0.  The  weighted  potential  is  most  easily 
obtained  by  integrating  pi  = p 94ii/3t. 

One  shot  from  the  experiments  of  Reference  7 was  analyzed.  The 
impact  velocity  of  the  round  was  1,354  ft/sec.  Inspection  of  the  high- 
speed motion  pictures  of  the  bullet  cavity  showed  that  the  bullet  had 
just  attained  a fully  tumbled  attitude  on  exiting  the  test  cell  . No 
permanent  deformation  of  the  exit  wall  was  observed,  so  the  assumption 
of  elastic  strain  in  the  analysis  is  justified. 

The  incident  pressure,  pi,  and  the  quantity,  p c uni,  are  plotted 
in  Figure  6 as  a function  of  time.  The  bullet  exits  the  test  cell  at 
1.6  milliseconds.  As  a consequence  of  the  diverging  flow  produced  by 
the  bullet,  pi  ^ p c Upi-  Problems  with  piston  theory  are  anticipated. 

The  exit  wall  deflection  was  then  calculated  by  the  variable  image 
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method  as  expressed  in  equation  15.  The  results  are  plotted  in  Figure  7, 
together  with  the  two  sets  of  experimental  points.  The  experimental 
error  is  estimated  at  ±0.25  cm.  The  variable  image  method  clearly  gives 
reasonable  agreement  with  experiment.  The  rise  time  of  the  experimental 
wall  deflection  is  much  faster  than  expected  from  theory.  This 
discrepancy  is  not  yet  understood. 

The  total  pressure  at  the  wall,  p,  and  the  source  strength,  k,  are 
plotted  as  a function  of  time  in  Figure  8.  The  incident  pressure,  p-j , 
is  included  for  reference.  During  the  first  part  of  the  panel  deflection, 
the  dominant  characteristic  of  the  plate  motion  is  inertia.  The  mass 
of  the  plate  is  so  low,  however,  that  the  plate  acts  very  much  like  a 
free  surface.  This  is  clear  from  Figure  8 where  p <<  p^ . The  source 
strength,  k,  responds  to  this  situation  by  dropping  very  rapidly  from 
its  initial  value  of  1 to  approximately  -0.9.  This  is  consistent 
with  the  previous  observation  that  the  variable  image  method  is  exact 
for  free  surfaces  when  k = -1 . 

As  the  wall  approaches  its  maximum  deflection,  the  bending  and 
membrane  stresses  begin  to  dominate  the  wall  motion.  The  inward 
acceleration  of  the  wall  is  resisted  by  the  fluid,  producing  a large 
net  pressure  at  the  wall.  The  image  strength,  k,  responds  to  this  by 
increasing  gradually  up  to  a value  of  3 or  4.  In  a sense,  the  wall 
might  be  called  "super  rigid"  at  this  point. 


FIGURE  6.  Weighted  Average  Incident  Pressure 
and  Normal  Velocity. 


The  variable  image  method  was  applied  here  assuming  a prescribed 
deflection  profile  of  the  exit  wall  and  a uniform  image  strength,  k. 
Equations  59  and  60  can  also  be  applied  when  using  more  sophisticated 
techniques  of  structural  analysis  by  allowing  k to  be  nonuniform  over 
the  wall . 
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Appciuiix  A 

DERIVATIVES  OF  POTENTIAL 
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Integral  in  00/dt  expression 
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PROGRAMMING  NOTE 


Tlie  remaining  integrals  are  difficult  to  evaluate  numerically.  For  example,  for  x term 
define 


|A(f)-B  t-tb(^) 


Find  an  expression  for  3I/9t  and  calculate  I as  a function  of  time  by  integrating  3I/3t. 


?ibnit/*s  Rule 
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Similarly,  for  the  cu  term 
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Appendix  B 

KLC  ()MMI:M)i;j)  liULLLT  CONSTANTS 


Conslanis 

.30  caliber 

.50  caliber 

1 2.7  mm 

14.5  mm 

HuJk't 

165.7  prams 

622-649  grains 

746  prams 

963  plains 

wcii;hl 

Hullol 

Piameter: 

Diameter; 

Diameter: 

Diameter; 

dinicnMons 

0.308  m. 

0.5I0->0.51 1 m. 

0.510-*  0.511  in. 

0.586  -♦0.587  in. 

Len>!lh:  1.39  in. 

Length;  2.31  in. 

Lenplh;  2.5  in. 

Length:  2.6  m. 

Pcnctratur 

80  -*81.2  prams 

375  grains 

wr4;ht 

448/449  prains 

655  grains 

Diameter: 

Diameter: 

Diameter: 

Diameter: 

dm\cnsuM\s  ' 

0.244  -►0.245  m. 

0.426  -►0.427  In. 

0.427  in. 

0.4889  in. 

Lenpth:  1.08  m. 

Length.  1 .8  in. 

Lenplh:  2.07  in. 

Length:  2.075  in. 

Jacket 

84.5  -►  85.7  prams 

247  -+  274  ptains 

297  -►298  prains 

293.21  grains 

wci):lil 

includes  incendiary 

includes  incendiary 

excludes  incendiary 

Jacket 

Same  as  bullet 

Same  as  bullet 

Same  as  bullet 

Same  as  bullet 

dimcnMi>ns 

dimensions 

dimensions 

dimensions 

dimensions 

BuIIct/jackcl 

projected 

area, 

untumBk'd 

0.0745  m^ 

0.2042  -►0.2050  m. 

0.2042  -►0.2050  in. 

0.2697  -0.2906  in. 

Bullcl/jjckct 

projected 

area, 

untumbled 

0.31732  m‘ 

0.7.3625  m^ 

1.03772  m2 

1.322  in2 

Penelrator 

projected 

area, 

untunibled 

0.0468  -►0.0471  m“ 

0.1425  -►0.1432  m^ 

0.1432  m2 

0.18773  in2 

Penelrator 

projected 

area, 

tumbled 

0.2601  m- 

0.59830  m2 

0.700214  m2 

1.05  in2 

Bullet 

projected 
area,  ittcrn- 
first 

0.0314  in^ 

0.1121  m2 

0.0855  m2 

0.12566  in2 

C()  of  ljullet, 

normal 

attitude 

0.10 

0.05 

0.05 

O.IO 

Ci)  of  bullet. 

tumbled 

attitude 

0.60 

0.30 

1 

0.30 

0.30 

Cj)  of  penelrator. 

normal 

attitude 

0.10 

0.05 

0.05 

0.10 

Cp  of  penetralor. 

tumbled 

attitude 

0.60 

0.30 

0.30 

0.30 

Cp  of  jaeket, 

tumbled 

attitude 

1.0 

1.0 

1.0 

1.0 

Cl)  of  bullet. 

Stern -first 
altiliide 

0.82 

0.82 

0.82 

0.82 

Degree  Obliquity,  0 Degree  Yaw. 


FIGURE  10.  Tumbling  Distance  Distribution;  0 Degree  Obliquity,  0 Degree  Yaw. 


TUMBLING  DISTANCE. 
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FIGURE  11  Tumbling  Distance  Distribution;  30  to  45  Degrees  Obliquity, 
0 Degree  Yaw. 
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